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THE  VALIDITY  OF  ASSUMPTIONS 
UNDERLYING  CURRENT  USES  OF 
LANCHESTER  ATTRITION  RATES 


STUDY 

GIST 


THE  REASON  FOR  PERFORMING  THE  STUDY  was  to  show  that  many  of  the  prevailing 
understandings  concerning  the  relationships  among  the  Lanchester,  Stochastic 
Lanchester,  and  the  General  Renewal  models  of  combat  are  erroneous  and  to 
collect,  organize,  and  reduce  to  common  notation  almost  all  known  tabled  results 
and  curves  of  particular  examples. 

THE  PRINCIPAL  RESULTS  of  this  study  were: 

(1)  All  Lanchester  model  and  Stochastic  Lanchester  model  mean  value 
equivalent  pairs  differ  for  all  times  except  at  possible  crossing  points. 
These  differences  may  be  very  large. 

(2)  At  least  for  the  Square  Law,  the  Lanchester  model  trajectories  are 
neither  a  universal  upper  or  lower  bound  of  the  Stochastic  Lanchester  model 
mean  value  trajectories. 

(3)  Even  near  time  zero,  the  Lanchester  and  Stochastic  Lanchester  model 
mean  value  trajectories  may  differ  considerably  (they  do  not  differ  materially 
for  the  Square  Law  and  the”  sequence  of  one-on-one  duels  version  of  the  Linear 
Law). 


(4)  For  the  Linear  Law,  the  Square  Law,  the  Mixed  Law  and  the  Square 

Law  with  continuous  reinforcements  there  is  a  Law  of  Large  Numbers  on  suitably 
transformed  spaces.  However,  for  untransformed  spaces  this  does  not  apply, 

for  it  can  be  shown  that  as  the  initial  force  sizes  tend  to  infinity  the 

differences  between  Lanchester  model  and  Stochastic  Lanchester  model  mean 
value  trajectories  tend  to  zero,  or  they  may  even  tend  to  a  constant  or 
infinity. 

(5)  Blackwell's  Theorum  does  not  imply  that  individual  combatants  with 

general  interfiring  times  tend  to  have  negative  exponentially  distributed 
interfiring  times.  This  is  even  more  strongly  the  case  for  terminating 

processes. 

(6)  The  Palm-Khintchine  Theorem  does  not  imply  that  superposing  a  large 
number  of  combatants  with  general  interfiring  times  will  yield  a  process  with 
negative  exponentially  distributed  interfiring  times.  This  can  only  be 
approximately  correct  for  large  numbers  and  for  very  large  interfiring  time 
means.  Again,  the  Theorem  is  only  valid  for  non-terminating  processes. 

(7)  Non homogeneous  Poisson  processes  do  not,  in  general,  approximate 
general  renewal  processes. 
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(8)  The  Stochastic  Lanchester  model  process  variances  are  generally 
quite  significant  and  can  be  important  for  large  force  sizes,  even  near  time 
zero.  In  addition,  general  renewal  model  process  variances  are  significantly 
different  that  Stochastic  Lanchester  model  process  variances. 

(9)  The  other  Lanchester  model  measures,  (a)  expected  number  of  survivors, 

(b)  expected  time  duration  of  the  battle,  and  (c)  probability  of  winning  are 
even  less  reliable  predictors  than  the  mean  value  trace. 

(10)  The  basic  assumptions  of  the  Stochastic  Lanchester  models,  as  well 
as  the  general  renewal  model,  cannot  hold  for  large  numbers  of  combatants. 

THE  MAIN  ASSUMPTIONS  were  that: 

(a)  All  pre-combat  decisions  have  been  made  and  the  battle  goes  forward 
until  terminated  by  the  action  itself  or  by  tactical  decisions. 

(b)  The  true  model  which  is  to  be  approximated  is  termed  a  General  Renewal 
model . 

(c)  Each  marksman  fires  until  he  is  killed  or  makes  a  kill. 

(d)  The  ammunition  supply  is  unlimited. 

(e)  All  fire  independently. 

THE  MAJOR  RESTRICTION  is  that  the  general  renewal  model  is  a  superposition 
of  many  terminating  renewal  processes. 

THE  SCOPE  OF  THE  STUDY  is  to  show  that  the  basic  assumptions  of  current  combat 
models  do  not  hold  for  large  numbers  of  combatants. 

THE  STUDY  OBJECTIVES  are  to  show  that  many  factors  cause  large  scale  battles 
to  be  a  number  of  simultaneous  and/or  sequential  smaller  scale  engagements 
and  that  use  of  current  deterministic  Lanchester  models  is  incorrect. 

THE  BASIC  APPROACH  Is  two  fold.  First  the  Lanchester  model  applications  are 
examined  to  show  fallacies.  Second  all  known  tabled  results  and  curves  of 
particular  examples  are  provided  to  support  the  theoretical  discussion. 

THE  STUDY  SPONSOR  is  TRAC-WSHR. 

THE  STUDY  PROPONENT  is  TRAC-HSMR. 


THE  ANALYSIS  AGENCY  is  TRAC-WSMR. • 
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NOTATION 


Where  multiple  definitions  are  given,  the  context  will  make  clear  which 
one  is  used. 

A  -  designation  of  one  side  In  the  combat,  or 

-  the  event  of  a  win  by  side  A 

Aj.(t)  -  rv,  the  number  on  the  A  side  at  time  t,  in  the  SL  or  GR  model  with 

ics  (kaQ,kbQ>  (when  k  =  l,  subscript  dropped),  jointly  distributed 

with  Bj.(t) 

a  -  value  of  rv,  A^(t) 

a^  -  the  number  on  side  A  at  the  time  the  A  side  loses  (breaks  and  runs) 
a^  -  the  initial  number  on  side  A  (at  time  zero) 

B  -  designation  of  one  side  in  the  combat,  or 

-  the  event  of  a  win  by  side  B 

Bj.(t)  -  rv,  the  number  on  the  R  side  at  time  t,  in  the  SL  or  GR  model  with 

ics  (kaQ.kb^)  (when  k  =  1,  subscript  dropped),  jointly  distributed 

with  Ak(t) 

b  -  value  of  rv,  Bk(t) 

bp  -  the  number  on  the  B  side  at  the  time  the  B  side  loses  (breaks  and  runs) 

bn  -  the  initial  number  on  side  B  (at  time  zero) 

cov(X,Y)  -  covariance  of  rvs  X,Y 

c  -  arbitrary  constant 

0  -  the  event  of  a  draw 

d  -  arbitrary  constant 

E [ X ]  -  expected  value  of  rv  X 

E [ X !  -  estimate  of  the  mean  of  rv  X 

e  -  arbitrary  constant 
c 

-  complementary  distribution  function  of  rv  X 
f  -  arbitrary  constant 
f'(t)  - 

f^(x)  -  pdf  of  rv  X 

GR  -  the  general  renewal  model 


-  v  - 


g^(t)  -  average  reinforcement.  raLe  on  the  i-^tl-  side,  t  =  A  or  B 
ic  -  initial  condition 
ift  -  interfiring  time 

iid  -  Independent,  identically  distributed  (used  with  rvs) 

k  -  positive  integer,  either  an  index,  or 

-  the  number  of  the  transition  step  in  sLate-space 

L  -  the  Manchester  model 

mA(t)  =  E[A(t)],  marginal  mean  value  function  of  A(  t ) 

mR(t)  =  E[B(l>],  marginal  mean  value  function  of  B(t) 

mo  -  portion  of  force  engaged  on  A  side  in  Springall  model 

N(t)  -  rv,  number  of  events  which  have  occurred  at  time  t  in  a 
GR  renewal  process 

n(t)  -  a  value  of  the  rv  N(t) 

n(t)  =  E(N(t)],  mean  value  function  of  Nft) 

iii  -  portion  of  force  engaged  on  B  side  in  Springall  model 

ned  -  negative  exponentially  distributed 

P(i)  -  probability  that  the  i  side  wins,  i  =  A  or  B 

P(i,t)  -  probability  that  the  i  side  has  won  bv  time  t,  i  =  A  or  B 

p  -  marksman's  kill  probability  or 

=  6  /(a+6 ) 

p(i,t)  -  pmf  of  rvs  A( t )  and  B(t)  respectively,  i  *  a,h 
p(a,t|A)  -  conditional  pmf  of  rv  A(l)  given  A  wins  by  time  l 

p(b,t|B)  -  conditional  pmf  of  rv  B( t )  given  B  wins  by  Lime  t 

Pj  -  the  constant  kill  probability  of  all  contestants  on  the  i—  side,  i=A  or  B 

pt(t)  -  time-dependent  kill  probability  of  all  contestants  on  the  i—  side, 

i  =  A  or  B 

P|c(a>f,»l)  ”  PfA^Ct)  *  a,  B^(t)  3  b,  ic(kaQ.kbfl)]  (for  k  =  1  subscript  is 
dropped) 

pdf  -  probability  density  function 


-vi- 


pmf  -  probability  mass  function 

q  =  1  -  p,  marksman's  failure  probability  or 
=■  a/(o+6 ) 

rhs  -  right  hand  side  of  an  equation 

r^  -  fixed  individual  firing  rate  for  side  i,  i  =  A  or  B 

r^(t)  -  time-dependent  individual  f ir ing  rate  for  side  i,  t  *  A  or  B 

r^(a,b,t)  -  i's  general  kill  rate  (other  side's  attrition  rate),  i  =  A  or  B 

r(y)  =  fz(y)/F^(y) ,  GR  model  individual  instantaneous  kill  rate 

rv  -  random  variable 

fo 

S.(t)  »u  1  ap( a , 0 , t ) 
a=l 

bo 

S-(t)  =  3  1"  bp(0 , b, t ) 
b=l 

SL  -  the  StochasLlc  Lanchester  model 
Tp  -  rv.  Lime  duration  of  combat 

Tl)[i  -  rv’  t^me  duration  of  comhai  given  a  win  by  i  =  A,B  or  a  draw,  i  =  D 
t  -  time 

t f  -  time  at  which  L  hatlles  terminate 
t|  -  time  at  which  iiil  renewal  event  occurs 
V|A(t)]  -  variance  of  the  rv  A(t) 

Wk(t)  -  rv,  superposed  interkilling  times 
X  -  rv,  marksman's  ift 

X.  -  rv ,  ift  of  each  member  of  the  side,  t  =  A  or  B 

Xk ( t )  -  rv,  a  transformed  version  of  A^(t),  k  =  1,2,... 

x  -  a  value  of  X^(t)  or 

-  a  value  of  X  or 

-  a  value  of  X| 

Xj  -  a  particular  value  of  x 

X2  ~  portion  of  initial  x  force  initially  engaged  in  Springall  model 
x(t)  -  the  number  on  side  A  at  time  t,  in  the  L  model 

xL( t )  -  solution  to  the  standard  L  Square  law  equations  with  ics  (aQ,bQ) 


-  v  i  i  - 


Y(t)  -  rv,  marksman's  backward  recurrence  time  aL  time  t 

Y^(t)  -  a  transformed  version  of  B^(t),  k  =  1,2,... 

y  -  a  value  of  Y^(t)  or 
-  a  value  of  Y(t) 

yj  -  a  particular  value  of  y 

y^  -  portion  of  initial  y  force  initially  engaged  in  Springall  model 
y(t)  -  the  number  on  side  R  at  time  t,  in  the  L  model 

y^(t)  -  solution  to  the  standard  L  .Square  Law  equations  with  ics  (aQ,bn) 
Z  -  rv,  embedded  killing  process  in  the  GR  marksman  model 


Aa  -  change  in  a 
AaU)  =  mA(t)  -  x(t) 

AR(t)=  mB(t)  -  y(L) 

A i  |  i  ( t )  _*  A  ^  ( t)  given  a  win  by  i ,  t  =  A,  B 

a  =  p./u,  A  A's  Individual  kill  rate  (attrition  coefficient  for  B  side) 

A  A  = 

B  =  pD/bD  A  B's  individual  kill  rate  (attrition  coefficient  for  A  side) 

D  D  — 

Y  -  non-comhat  attrition  coefficient,  B  side  or 

-  Springall  attrition  coefficient 

A(x-a)  -  Dirac  Delta  function  =  00  at  x  =  a 

=  0,  elsewhere 

<5  -  a  constant,  or 

-  Springall  attrition  coefficient 

c  -  an  arbitrary  positive  constant 
5  -  mean  of  rv  X 

d  -  y  -  y{ 

Pq  -  an  arbitrary  positive  constant  <  1 

9  -  the  probability  that  an  A  side  target  is  not  acquired  by  a  B  side  firer 
A  -  the  probability  that  a  B  side  target  is  not  acquired  by  an  A  side  firer 

V  -  c/p  -  E[z] 


-  v  i  i  i  - 


II  t> 


U  -  mean  value  of  the  1ft  on  side  i  for  the  SL  or  GR  model,  1  *  A  or  B 
moment  about  the  origin,  k  =  1,2,... 

5  =  x  -  Xj 

5^  -  an  arbitrary  positive  constant  <  1 
P  -  non-combat  attrition  coefficient,  A  side 

a  -  standard  deviation  of  the  losses  on  the  A  side  in  the  SL  process  at  time 
of  the  k1-^  event 


- +  -  converges  in  distribution  to 

p 

- ♦  _  converges  in  probability  to 

*  -  cartesian  product 

defined  as  equal  to 


-  i  x- 
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PART  ONE  -  THEORY 


I.  INTRODUCTION 

Great  emphasis  Is  placed  In  this  document  on  assumptions,  results,  and 
concepts  as  they  relate  to  small-to-moderate  size  battles.  The  reason  for 
this  Is  that  the  basic  assumptions  of  current  combat  models  do  not  hold  for 
large  numbers  of  combatants.  Terrain  compartmentallzation,  weapon  ranges, 
terrain  obstacles,  weather  and  many  other  factors  cause  large  scale  battles  to 
be  a  number  of  simultaneous  and/or  sequential  smaller  scale  engagements. 

Much  of  the  current  effort  in  weapons  systems  analysis  or  in  combat 
analysis  relies  on  utilization  of  various  deterministic  Lanchester  (L)  models 
of  combat.  It  Is  important  to  emphasize  at  once  that  the  models  discussed 
here  are  only  concerned  with  the  progress  of  the  fire-fight  once  it  has 
begun.  All  pre-combat  decisions  have  been  made  and  the  battle  goes  forward 
until  terminated  by  the  action  Itself  or  by  tactical  decisions  (e.g.  withdraw, 
surrender,  etc.).  This  means  that  optimization  techniques  which  influence 
pre-combat  tactical  decisions  or  pre-determined  decisions  to  terminate  are 
excluded.  Optimization  employs  Differential  Game  Theory,  Dynamic  Programming, 
Game  Theory,  Control  Theory,  certain  mathematical  programming  procedures  or 
possibly  other  optimization  techniques  and  requires  good  fire-fight  models  to 
operate  on. 

The  most  critical  element  in  fire-fight  analyses  is  the  assignment  of  the 
attrition  coefficients.  For  gxample,  the  Manchester  Square  Law  differential 
equations  (see  reference  {101  )  governing  the  solution  are 

x'(t)  =  -fiv,  y ' (t)  =  -nx, 

where  x  and  y  are  taken  to  be  the  average  numbers  remaining  on  the  two  sides 
(A  and  B)  respectively  and  where  a  and  8  are  the  B  and  A  sides'  attrition 
coefficients  (A  and  B  individual  kill  rates)  respectively.  It  is  customary  to 
use  p^/u .  for  o,  where  p^  is  the  individual  single  round  kill  probability  on 
the  A  side  and  is  the  individual  mean  Interfiring  time  on  the  A  side.  All 
individuals  on  each  side  are  assumed  to  be  identical.  \  similar  set  of 
assumptions  and  definitions  go  with  the  8  side. 

A.  The  Assumptions  Implicit  In  Using  L^ncheater's  Equations 

The  analyses  In  this  paper  are  based  on  the  following  additional 
assumptions,  which  are  explicated  for  the  Square  Law  for  the  sake  of  specificity. 
The  details  of  other  models  will  vary,  hut  the  analysts  Is  similar  and  the 
attrition  coefficients  are  defined  in  a  manner  appropriate  to  the  model  being 
considered.  These  assumptions  are  the  basis  for  using  L  models  mentioned  above. 

(1)  The  true  model  which  is  to  be  approximated,  is  termed  a  General  Renewal 
(GR)  model,  and  is  Illustrated  in  Figure  1  below.  The  characteristics  of  this 
model  are: 

(a)  There  are  initially,  Hq  on  the  A  side  and  b^  on  the  B  side. 

(b)  Every  member  of  the  A  side  picks  a  B  opponent  at  random  (all  are 


*  Numbers  in  brackets  [  ]  refer  to  the  list  of  references  starting  on  p.  40. 
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visible  and  in  range)  and  fires  with  a  fixed  kill  probability  on  every 
round  fired  with  a  general  interfiring  time  (ift),  which  is  a  random 
variable(rv)  XA  with  mean  u..  (Because  a  general  rv  has  memory  some 
particular  distribution  of  fire  is  required  for  the  GR  model  to  be  formulated 
mathematically.  Different  distributions  of  fire  will,  in  general,  give 
different  results.  However,  no  matter  what  distribution  of  fire  is  assumed  in 
the  GR  model,  as  long  as  all  combatants  have  a  target  (even  if  all  are  firing 
at  the  same  target  for  example)  the  GR  model  reduces  to  the  SL  model  when  XA 
and  Xg  are  assumed  negative  exponentially  distributed  (ned).  Thus  we  see  that 
for  L  and  SL  the  restriction  that  all  opponents  are  visible  and  in  range  can 
be  reduced  to  simply  that  every  firer  can  see  and  possibly  hit  at  least  one 
target  at  all  times.)  For  the  GR  model,  it  must  be  specified  not  only  how 
many  opponents  are  targets  but  how  the  fire  is  distributed  over  targets.  This 
makes  for  possibly  many  different  GR  models  for  every  SL  and  corresponding  L 
model.  To  our  knowledge  this  has  not  been  pointed  out  before. 

(c)  Each  marksman  fires  until  he  is  killed  or  makes  a  kill  (at  which 
point  he  Immediately  shifts  to  a  new  target  picked  at  random  and  resumes 
firing).  The  ammunition  supply  is  unlimited. 

(d)  All  fire  independently. 

(e)  Similar  assumptions  apply  to  the  B  side. 

(f)  The  battle  continues  until  one  side  is  annihilated. 


A  side  8  side 


1 

^-0 

• 

• 

• 

• 

• 

0^ 


Figure  l.  The  GR  Model 

The  derived  quantities  of  Interest  are,  A(l)  and  B(t),  rvs  which  are  the 
remaining  numbers  on  each  side  with  expected  values  F.fA(t))  =  m^(t)  and 
EfB(t)]  =  mp(t).  We  note  that  these  are  marginal  rvs  and  means  from  a  joint 
distribution  on  some  state-space. 

This  situation  is  a  superposition  of  many  terminating  (non-classical) 
renewal  processes  (i.e.  each  marksman's  firing  process  Is  a  renewal  process 
which  terminates  but,  in  general,  not  at  one  of  his  firing  epochs).  The 
reader  is  reminded  that  an  ordinary  general  renewal  process  is  one  in  which 
the  times  between  events  are  Independent,  identically  distributed  (lid)  rvs 
and  that  never  terminates. 
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The  only  solutions  to  this  model  have  been  for  the  one-on-one  duej^  (see 
reference  [1])  and  for  the  two-on-one  duel  (fiafarian  and  Ancker  (1984)  ). 

(2)  It  is  assumed  that,  if  aQ  and  are  verv  large,  mA(t)  and  mR(t) 
are  well  approximated  by  the  mean  value  solutions  to  the  “Stochastic 
Lanchester"  (SL)  equations  which  were  first  given  a  thorough  treatment  by  Snow 
(1948).  The  SL  model  is  exactly  the  GR  model  except  that  ifts  are  assumed  to 
be  negative  exponentially  distributed  (ned).  The  ned  assumption  with  its  “no 
memory"  property  greatly  simplifies  the  analysis  hut  does  not  make  it  trivial. 

(3)  It  is  assumed  that  for  large  a^  and  b^,  the  L  model  solution  is  a 
good  approximation  to  mA(t)  and  mg(t)  ( the  mean  value  time  traces)  of  the  SL 
model.  Lanchester  (see  reference  [10])  firsL  proposed  his  model  in  1916  and 
recognized  that  since  it  was  deterministic  and  continuous  in  x,y  and  that  the 
real  process  was  stochastic  and  discrete,  L  should  be  considered  an 
approximation  of  the  stochastic  process  mean  value. 

(4)  Finally,  it  is  assumed  that  for  large  a^  and  Hq  the  SL  and  GR  models 
have  negligible  variance  for  much  of  the  battle. 

It  is  noted,  parenthetically,  that  we  consider  all  these  models,  (L,  SL 
and  GR)  to  be  logically  subsumed  under  the  title  "Lanchester"  models.  Thus  SL 
and  GR  refine  the  original  concepts  of  Lanchester  in  the  direction  of  reality 
at  the  obvious  cost  of  greater  complexity. 

It  is  this  set  of  assumptions,  (I)  through  (4)  above,  which  are 
thoroughly  examined,  particularly  the  basic  premises  which  purport  to  justifv 
them.  The  main  text  (Part  One)  of  this  paper  will  examine  all  known 
theoretical  results  on  this  subject  and  in  Part  Two  we  collect,  organize  and 
reduce  to  common  notation  almost  all  known  tabled  results  and  curves  of 
particular  examples.  These  support  the  theoretical  discussion. 

The  theoretical  discussion  will  proceed  in  reverse  order  from  assumption 
(3)  backwards  to  (2)  then  take  up  (4)  and  close  bv  considering  other  measures 
of  effectiveness.  First,  several  preliminary  matters  will  he  examined. 

B.  The  Lanchester  Square  Law 

It  will  he  helpful  to  look  at  the  solution  to  the  Lanchester  Square  Law 
differential  equations  given  above.  These  are  the  parameteric  equations 

x ( t )  =  a^  cosh  /<iB  t  -  /(’./a  b(1  sinh  * -v-'  t  , 

y(t)  =  cosh  if  all  t  -  /a/R  sinh  /nR  I  , 

and  the  phase-space  equation  (which  is  arrived  at  bv  dividing  the  second 
differential  equation  above  into  the  first  am)  solving), 

y(L)  -  {(a/8)  (x2(t)  -  a2!  +  b(2}  1/2 
all  of  which  are  only  valid  in  the  interval  (O.tp)  where. 


References  given  by  name(s)  followed  by  a  date  are  listed  in  the  Annotated 
Bibliography  beginning  on  p.  42.  if  the  date  is  replaced  by  an  asterisk  (*), 
the  material  Is  new  In  this  work. 
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survivors 


tf  -  (1/VafT)  tanh  ^^aTs  ybQ),  if  ''BA*  b0/aQ  >  1, 

-  (I/V06)  canh-1(/MT  b^),  if  ^  b0/aQ  <  1, 

-  if  bQ/a0  -  1, 


2  2  1/2 

B  wins,  with  (bg  -  a^x/B) 

2  2,  L/2 

A  wins,  with  (a^  -  bj/a)  survivors 
which  is  a  draw  with  no  survivors. 


At  time  t^,  the  combat  has  finished  and  one  (or  both)  side(s)  is  annihilated. 
Beyond  that  time,  the  only  interpretation  which  can  be  useful  is  to  assume  the 
states  remain  frozen  at  their  t^  values.  This  is  illustrated  in  Figure  2  below. 


Figure  2.  Typical  Solutions  to  the  Lanchester  Square  Law  Equations 
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The  important  point  to  observe  here  is  that  at  tf  the  equations  are  no 
longer  valid  as  they  show  an  upturn  on  one  side  (impossible  as  there  are  no 
reinforcements)  and  at  exactly  the  same  time  the  other  side  starts  to  go 
negative  (again  an  obvious  impossibility).  The  only  exception  to  this  is  when 
parity  exists  (i.e.  ag  /a  ■  bg  /B  );  then  t^  is  infinite  and  the  situation  is 
as  depicted  in  Figure  2(c). 

C.  Properties  of  the  Equivalent  Stochastic  Solutions 

It  is  useful  to  recall  some  properties  of  the  joint  probability  mass 
function  (ptnf)  of  A(t)  and  B(t)  (these  apply  equally  well  to  the  SL  or  GR 
models). 

On  the  time-evolution  joint  probability  state-space  diagram  (see  Figure 
3),  the  only  possible  states  are  at  the  lattice  points  (the  intersections  of 
the  coordinate  lines  shown).  At  time  t,  the  probability  of  being  at  a  general 
lattice  point  (a,b)  is  designated  by  p(a,b,t).  The  states  along  the  a-axis 
and  the  b-axis  are  absorption  states  (all  opponents  have  been  eliminated). 

All  other  states  are  transient.  The  state  (0,0)  can  never  be  occupied.  At 
time  zero,  the  process  is  at  (ag,bg)  with  probability  one.  For 


Figure  3.  The  Time-Evolution  Joint  Probability  State-Space 

all  t  >  0,  all  allowable  states  have  non-zero  probability  of  occurring.  The 
point  (m^( t ) ,mg( t ))  is  where  the  mean  value  function  pierces  the  state-space 
plane.  In  the  limit  as  t  ♦  •»,  all  transient  states  have  zero  probability  of 
occurring  and  all  the  probability  mass  is  located  on  the  a  and  b  axes.  This 
Is  still  a  joint  pmf  but,  of  course,  all  probabilities  are  just  on  the  two 
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axes  as  shown.  The  probability  the  A  side  has  won  by  time  t  is  just  the  sum 
of  all  the  probabilities  on  the  a-axis,  and  similarly  for  the  B  side.  The 
total  win  probabilities  are  these  functions  evaluated  at  t  *  «  . 

There  seems  to  be  some  confusion  in  the  literature  about  the  joint  prnfs 
shown,  the  marginal  pmfs  and  the  conditional  absorption  prafs.  In  Figure  3  (at 
any  time,  t)  the  marginal  pmf  for  A  is  simply  all  the  probability  mass 
projected  onto  the  a-axis  and  similarly  for  B.  Given  that  A  has  won,  the 
conditional  pmf  is  found  by  dividing  each  mass  point  on  the  a-axis  by  the  sum 
of  the  mass  points  on  the  a-axis  and,  similarly  for  B.  Particular  confusion 
appears  to  occur  in  connection  with  moments  of  the  three  possible  prafs  (joint, 
marginal  and  conditional)  especially  in  the  terminal  state  (t  =  ®)  and  it  Is 
important  to  keep  clearly  in  mind  which  is  intended. 

On  the  time-independent  joint  probability  state-space  (see  Figure  4)  the 
diagonals  labelled  k  *  1,2,...,  (ag+bg-l)  contain  the  possible  states  after  the 
first,  second,  third,  etc.  events  have  occurred.  The  probability  masses  located 
on  each  diagonal  (k  <  min  (ag,bg))  form  a  proper  pmf  and  add  to  one.  For  k  > 
min  (ag,bg),  all  absorption  probabilities  associated  with  smaller  ks  must  be 
included  with  those  on  the  diagonal  to  form  a  proper  pmf  (see  e.g.  the  broken 
line  for  k  =  ag  in  Figure  4.)  A  random  walk  may  be  visualized  on  this  space  with 
each  step  either  to  the  left  or  downward  and  terminating  on  the  a  or  b-axis. 


Figure  4.  The  Time-Independent  Joint  Probability  State-Space 

It  is  important  to  note  here  that  the  concept  of  "breakpoint"  (i.e., 
defeat  before  annihilation)  complicates  the  mathematics  but  does  not  change 
the  conclusions  of  the  following  analyses  in  any  material  way.  Breakpoints 
create  absorbing  barriers  at  specific  positive  values  of  a  and  b  (ap  less  than 
ag  and  bp  less  than  bg).  This  means  that  absorption  probabilities  become 
significantly  large  earlier  in  the  battle  and  their  role  In  the  differences 
between  L  and  SL  measures  of  battle  outcomes  Is  effective  earlier. 

We  now  proceed  to  examine  important  theoretical  points  about  the 
assumptions  outlined  above. 


II.  THE  ERROC  IH  COgSIDERIMG  LAHCHESTER'S  SOLUTIONS  AS  AH  APPROXIMATION  TO 
STOCHASTIC  LAHCHESTER  MEAN  VALUE  FUNCTIONS 


Before  proceeding  with  details,  it  is  noted  that  essentially  the  question 
here  is  -  how  well  does  the  L  joint  trace  (e.g..  Figure  2(d))  approximate  the 
corresponding  SL  mean  value  trace  of  Figure  3  (or  equivalently  their  marginal 
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or  phase-space  projections)? 

In  an  Important  paper  not  readily  accessible  and  thus  largely  unnoticed 
Hardeck  and  Hllden  (1967,  equations  (10),  p.  5)  have  shown,  when  their  result 
Is  expressed  In  terms  of  the  number  of  survivors,  that 


mA(t) 


mB(t) 


a0  bo* 

l  l  r„(a,b,t)  p( a, b, t )  =  -£(r_(A(t),B(t),t) J 

a^a^+l  b*b^+l 

a0  b0 

l  l  r, (a,b,t )  p( a , b ,  t )  =  -E  f  r  (A(t ) , B( t ) , t )  ] 

a=af-H  b=bf+l 


O) 


for  a  general  SL  process  (with  breakpoints)  and  where  kill  rates  (r^,r^)  also 
depend  on  a  and  8  even  though  this  dependency  is  not  explicitly  shown.  Later, 
for  the  same  general  SL  process  with  kill  rates  not  a  function  of  time  and 
where  the  combat  goes  to  annihilation,  Clark  (1969  equations  (73),  (74)  and 
(75),  p.  79)  has  derived  these  same  expressions  except,  of  course,  that  the 
kill  rates  were  not  dependent  on  time. 


Although  neither  Hardeck  and  Hi Iden  nor  Clark  sav  so,  these  expressions 
conclusively  prove  that  compared  to  any  equivalent  L  formulation,  the  SL  form 
will  be  different. 


This  comes  about  because  the  right  hand  side  (rhs)  of  equation  (1)  cannot 

be  made  to  look  like  the  rhs  of  equivalent  L  equations.  There  are  two  reasons 

this  cannot  he  done.  They  are: 

(1)  There  are  no  absorption  probabilities  utilized  in  the  SL  formulation  for 

computing  the  expected  kill  rates  because  kill  rates  are,  by  definition,  zero 

as  soon  as  one  side  is  annihilated  (or  reaches  Its  breakpoint). 

(2)  A(t)  and  B( t)  are  correlated. 

in  any  possible  L  model  one  or  the  other  or  both  of  these  facts  will  create  a 
difference  between  the  L  differential  equations  and  the  equivalent  SL  mean 
value  differential  equations. 

This  will  be  illustrated  below  by  all  particular  examples  in  the 
literature  known  to  the  authors. 

A.  The  Square  Law 

1.  Time-Independent  Kill  Rates 

The  first  analytic  attention  to  this  problem  was  given  hy  Snow  (1948). 
However,  he  first  gives  a  more  general  L  problem  (see  equations  (3)  and  (25) 
in  Snow);  namely 
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and 


x'(t)  -  -Px  -  3y  +  j?AC t > ,  | 
y'(t)  -  -Yy  -  ax  +  gg(t),  ) 

then  derives  (see  p.  24,  Snow)  the  more  general  SL  equivalent 


m^(t)  *  -PmA  "  +  +  3  1  bp(0,b,t), 

b=l 


m^(t)  =  -YmR  -  ct^  +  gB(t)  +  a  l  ap(a,0,t)  , 

a=l 


(2) 


where  p  and  Y  are  non-combat  attrition  coefficients  and  g^(t)  and  gB  ( t )  are 
reinforcement  functions. 


It  should  be  observed,  in  passing,  that  Tompkins  (1953,  p.  37,  equation 
3.2)  rederived  these  results  (with  gA(t)  =  gR(t)  5  r>),  apparently  being 

unaware  of  Snow's  priority. 

The  Square  Law  equations  are  easily  derived  bv  setting  p  =  y  =  g^Ct)  - 
g  (t)  =  0  in  equations  (2)  and  (3).  Mote  that  the  L  differential  equations 
differ  from  the  SL  differential  equations  by  the  Lerms  containing  absorption 
probabilities.  The  SL  Square  Law  is  also  easily  derived  from  equation  (1)  Sv 
setting  r^(a , b, t )  ■  rA(a,b)  =  aa  and  rR(a,b,t)  =  rR(a,b)  =  8 b. 

We  note  from  the  general  differential-difference  state  equations  (Clark 
1969,  equations  (51)  through  (57)  pp.  69,70)  with  Square  Law  rates  inserted 
that  the  state  probabilities  have  the  following  properties  (see  Figure  5): 

(1)  The  absorption  probabilities  (p(a, 0,0,  p(h,0,t))  shown  in  Figure  5(i) 
are  monotonical ly  increasing  functions  with  value  zero  and  slope  zero  at  t  =  0 
and  asymptotes  p(a,0,°°)  and  p(0,b,°°).  The  curvature  is  initially  positive 
and  changes  to  negative. 

(2)  The  Initial  state  probability,  n( an , bn , L  ) ,  is  a  mono ton tea 1 1 v  decreasing 
function  with  value  one  and  with  slope  -(a^a  +  h^0  )  at  t  =  0  and  is 
asymptotic  to  zero  at  infinity.  The  curvature  is  alwavs  positive  (see  Figure 
5(b)). 

(3)  All  other  state  probabilities  have  values  zero  and  slope  zero  at  time 
zero  and  are  asymptotic  to  zero  at  infinity.  The  curvature  is  initially 
positive,  changes  to  negative  then  goes  back  to  positive  (see  Figure  5(c)). 

These  facts  along  with  the  equations  for  the  mean  value  functions  for  the 
Square  Law  show  that  m^(t)  and  mB(t)  have  values  iQ,bQ  and  slopes  -Bb^  and 
-aaQ,  at  t  *  0  and  are  monotonically  decreasing  to  positive  asymptotes  with 
slopes  zero  at  infinity.  They  also  have  positive  curvature  everywhere  tending 
to  zero  at  infinity. 

The  Square  Law  L  equations  have  the  same  initial  values  and  slopes  as  the 
SL  mean  value  functions  and  have  positive  curvature  (hut  different  from  SL) 
for  positive  t  <  t^.  It  should  he  noted  that  in  comparing  all  L  and  SL  pairs 
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1 


»<o,0,»l  or  p(0,b,«) 


Figure  5.  Characteristics  of  SI.  Square  Uw  State  Probabilities 

of  equations,  x  Is  always  identified  with  and  v  with  m,j.  This  Is  Lhe  onlv 
rational  interpretation  of  x  and  y  and  has  always  been  so  recognized  by  most 
analysts  including  Lanchester  himself  as  mentioned  earlier.  Also  iL  should  be 
observed  that  mA(t),  m^(t),  x(t)  and  v(t)  will  be  written  as  m^,  mR,  x,  and  v 
on  the  rhs  of.  all  equations  for  brevity,  always  keeping  in  mind  that  they  are 
functions  of  time. 

With  these  conventions  in  mind.  It  is  clear  that  equations  (2)  and  (3) 
differ  only  in  that  the  rhs  of  the  SL  equations  contain  terms  proportional  to 
the  average  number  of  survivors.  Clark.  (I9f>9)  has  dubbed  these  the  "bias" 
terms  of  the  SL  equations.  This  terminology  is  dropped  In  this  paper  because 
bias  has  a  well  known,  rigorous  meaning  In  statistics  which  does  not  ippiv 
here  as  this  paper  deals  with  purely  probabilistic  models  with  no  question  of 
statistical  sampling  Involved.  The  focus  here  is  on  the  dif Ferenee  in  the 
solutions  to  the  differential  equations  (l.e.  the  difference  between  the  i. 
functions  and  SL  mean  value  Functions). 

The  principal  point  to  note  here  Is  thaL  for  t  >  0 

l  >  p(a,0,t)  >  0  ,  0  <  a  <  nQ  > 
l  >  p(0,b,t)  >0  ,  0  <  b  <  b0  , 

and  these  probabilities  are  monotonical ly  increasing  (their  derivatives  are 
positive  for  all  t).  Thus,  although  mA(t)  and  x(t),  and  mR(t)  and  y(t)  start 
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at  the  same  point  at  t  *  0  (their  Initial  conditions  require  raA(0)  **  x(0)  =  ag 
and  mg(0)  *  y(0)  =•  bg)  in  general,  the  L  and  SL  mean  value  trajectories  will 
differ  for  all  t  (the  only  exception  is  that  they  may  cross,  and  therefore  be 
equal,  at  certain  specific  times). 

2.  Tiae-De pendent  Kill  Rates 

Hardeck  and  HiLden  (1967,  p.  8)  have  shown  that  for  time-dependent  kill 
probabilities  pA(t)  and  pB(t)  and  firing  rates  rA(t)  and  rg(t)  that 


0 

mA(t)  =  -pB(t)rB(t)  [mB  “  l  b  p  (  0 ,  b ,  t )  j  , 


n>g( t)  =  -pA(t)rA(t)[mA  -  I  a-  p(a,0,t)]  . 

a=  1 

The  corresponding  L  equations  again  only  differ  by  the  terms  containing  the 
absorption  probabilities,  and  the  comments  in  section  A.  1.  above  apply. 


This  is  fundamentally  different  than  the  other  processes  that  are 
examined  here  in  that  the  kiLLing  process  embedded  in  the  firing  process  is 
nonhomogeneous  Poisson  and  is  therefore  not  renewal  (i.e.,  it  is 
not  lid). 


B.  The  Square  Law  with  Breakpoints 

Craig  (1975  equations  (101),  and  (102)  pp.  160  and  161)  has  shown 


J0 


*0 


AU)  =  -0mB  +B  {  [  bp(af,b,t)  +  bf  l  p(a,bf,t)}, 


b=bf+l 


a"Vl 


(4) 


rn^(t)  =  -amA  +  a  {  £  ap(a,bf,t)  +  nf  l  p(af,b,t)f. 


a=  a  J.  +  L 


b=bf +1 


The  corresponding  L  equations  are  the  same  as  the  annihilation  case  except 
that  ap  and  bp  occur  in  the  boundary  conditions  in  an  obvious  way. 


Although  the  sums  on  the  rhs  of  equations  (4)  (which  produce  the  L-SL 
differences)  are  more  complicated  than  for  equations  O),  the  conclusions  are 
not  altered  materially.  Again  we  get  the  simple  Square  Law  by  letting  ap  =  bp 
=  0. 


C.  The  Linear  Law 


For  the  well  known  L  Linear  Law  given  by, 

x'(t)  -  -Bxy,  y'(t)  =*  -axy,  (5) 


-10- 


y  »  (a/6)x  +  (Bbg  ~  aao^®‘ 

AIL  these  solutions  are  valid  In  the  interval  ( 0 ,  t ^ )  where  tp  =  00  .  Tf 

aaQ  >  Sbpj,  A  wins,  x(®)  =  (aaf)  “  Bb())/a,  /(“)  =  0, 

aa0  <  8b^,  B  wins,  x(°°)  =  0,  y(°°)  =  (S -  aa^/S, 

aa^  =  Bb^,  draw  with  no  survivors,  xC00)  =  y(°°)  =  0. 

dark  (1969  equations  (80)  and  (81)  on  p.  81)  has  shown  that  the  equivalent 
SL  differencial  equations  are 

m'(t)  =  -BE [ A( t )B( L ) ]  =  -8m  m  -  8  cov [ A( L )B(  L ) I  , 

A  AH 

nrig(t)  =  -aE  [  A(  t)  B(  t )  ]  =  -amAmR  “  a  cov  ( A(  t )  B(  t )  ] 


It  should  be  noted  here  that  the  phase-space  solution  for  equations  (8) 
derived  by  dividing  the  first  equation  by  the  second  equation  and  solving  are 
exact ly  the  same  as  obtained  by  the  same  process  in  equations  (6).  Other 
versions  of  the  Linear  Law  will  give  the  same  phase-space  equations,  but 
different  time  traces.  In  general,  replace  xy  on  the  rhs  of  equations  (5)  bv 
any  general  function  g(x,y)  and  a  type  of  Linear  Law  will  result.  The  fact 
that  the  L  and  SL  time-independent  phase-space  equations  are  exactly  the  same 
(a  straight  line  starting  at  x(0)  =  3q,  v(0)  =  b^  and  m^(0)  =  ap,,  mR(0)  - 
and  with  the  same  slopes)  does  not  mean  the  lines  are  Identical,  The  L  lines 
terminate  either  on  one  of  the  axes  or  the  origin,  the  SL  lines  terminate 
at  (m^(°“ )  ,m_(® ) )  which  are  always  positive  (see  Figure  (6a)).  'Inless 
otherwise  stated  when  "Linear  Law"  Is  used  it  will  mean  the  versions  given  Ln 
equations  (5)  and  (6). 

Clark  did  not  actually  show  the  second  version  of  the  rhs  of  equations 
(6),  but  they  are  obvious  from  basic  probability  theory.  Again  the  L  and  the 
SL  mean  value  differential  equations  differ,  this  time  by  the  covariance 
terms.  Thus  the  SL  version  could  only  be  the  same  as  the  L  version  If  A(t) 
and  B(t)  were  Independent  or  uncorrelated.  Since  all  the  p(a,h,t)  functions 
with  a,b  *  0  are  known  (Clark  (1969),  p.  102  equation  (106))  and 
are  not  of  the  form  p(a,t)  p(b,t)  (l.e.,  product  of  the  marginal 
probabilities),  then  A(t)  and  B(t)  are  not  Independent  and  therefore  A(t)  and 
B(t)  may,  but,  In  general,  will  not  have  zero  covariance.  Thus  although  the 
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OL  Equation  Terminal  Point 
*  SL  Equation  Terminal  Point,  (mA(co), mB(ao)) 


Figure  6.  Comparison  of  L  and  SL  Solutions  for  the  Linear  Law 

phase-space  equations  are  eolllnear,  the  mean  value  time-traces  for  L 
dLffer  (see  Figure  6b). 

Equation  (6)  may  easily  be  derived  from  equations  (l)  by  letting 
rA(a»b,t)  *  rA(a,b)  »  aab  and  r^Ca.b.t)  =  r^Ca.b)  =  Bab. 

D.  Two  Special  Models 

(l)  Clark  (1969,  p.  151)  has  Investigated  a  special  case  where 
acquisition  probabilities  are  involved.  The  L  equations  are 


and  SL 
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x'(t)  *  -8y(l-0x)  =  -8y  +  8y9X 
y'(t)  =  -ax(l-A^)  =  -ax  +  axA-V 

and  the  equivalent  SL  equations  are 

m'(t)  -  -6E[B(t)(l-9)A(t)]  =  -8m  +  Bm„S 

A  B  B 

+  8  co 

m^(t)  »  -aE  [A(  t )(  l -1  ^  ]  =  -am^  +  l*mRE[; 

+  o  cov 

where  0  <  9  ,A  <  1  are  target  nonacquisition  probabilities  for  A  and  3 
respectively. 

Again,  there  is  a  difference  in  the  two  versions  which  is  even  more 
A  mA 

pronounced  as  E[9  ]  cannot  yield  9  .  We  note  that  equation  (8)  may  be 

obtained  from  equations  (l)  by  letting  r,  =  aafl-A^)  '  rn  =  Bb(l-9a). 

A  IS 

(2)  Springail  (1968)  has  thoroughly  investigated  a  rather  complicated 
model  whose  L  formulation  is  given  bv, 

x'(t)  =  -6xy  -  <5v, 

y’(t)  =  -axy  -  yv, 

where  5  and  Y  are  additional  fixed  attrition  coefficients.  The  initial 
numbers  engaged  are  Xt  and  y2 ,  which  are  fractions  of  the  initially  ava i lable 
forces  (aq,  b q  respectively).  The  remaining  forces  are  In  reserve  and  are 
deployed  one  by  one  as  the  initally  engaged  forces  are  reduced  in  such  a 
manner  as  to  keep  the  engaged  forces  at  levels  x0  and  v,  respect ive Iv ,  until 
all  reserves  are  committed  and  then  the  battle  proceeds  either  to  annihilation 
or  to  a  specified  breakpoint. 

The  corresponding  SL  mean  value  functions  ire  too  complicated  to  be 
reproduced  here  but  examples  of  the  outcome  of  such  battles  are  given  In 
Part  Two  and  show  substantial  L  function  -  SL  mean  value  function  differences. 

E.  An  Analysts  of  the  L-SL  Mean  Value  Function  Difference 

Define  the  L-SL  mean  value  function  differences  to  be 


(7) 


•  ) 

[9A(t)] 

v[B(t)0A(t)]  , 
[A(t)AR(°l  , 


AA(t) 

AB(i) 


x  , 


mB  -  y  . 


I 

I 


(9) 
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1.  The  Square  Law 


Craig  (1975,  p.  68)  has  shown  that  for  the  SL  Square  Law 


A^( t )  =  /q  {Sb(t)  cosh  /aS  (t-t) 
Ag( t)  =  /q  {Sa(t)  cosh  /a8  (t-t) 
where  0  <  t  <  and 

an 

SA(t)  =  a  l  ap(a,0, t ) , 
a=l 


S_(t)  -  8  l  bp( 0 , b , t ) . 
8  b=l 


i /8/a  Sa(t)  sinh  /a6(t-T)}dr 
/a/8  S  (t)  sinh  /aB  (t-i)}dx 


(10) 


SA(t)  and  S^(t)  are  the  terms  in  the  SL  dif ferent ia 1  — d i f f erence  equations  for 
the  ordinary  Square  Law  formulation  causing  L-SL  differences.  StricLlv 
speaking,  Craig  uses  the  breakpoint  model  for  equations  (9)  and  (10)  but  this 
includes  the  annihilation  model  we  are  considering.  Craig  investigated  (10) 
by  using  various  particular  values  of  an,  b^,  a,  8.  We  shall  go  a  bit  further 
here  to  draw  some  more  general  conclusions. 


Equations  (10)  mav  easllv  be  rewritten  as 


AA(t)  =  {/J  [Sb(t)  -  /eTo-  S a ( t  )  |  exp f  / a8  (t-r)ldT  \ 

+  /or-VT>  +  ,/S//a  sA(t)|  exp  1/(16  ( t-r  1  |  dr  } /2,  I 

)( I  I  ) 

Wow,  by  the  Second  Mean  Value  Theorem  equation  (II)  can  be  written  as 


A A ( t )  -  { l-exp(-  /«8t )} {  1 1 (c, t )  +  I2(d,t)}  /2/a8  ,  \ 

AR(t)  =  {  l-exp(-  /aBt )}  {  -I .  (c  ,  t )  +  I2(d,t)}/28,  whore  I 

_  _  I  (12) 

I[(c,t)  =  [SR(c)  -  /8/a  SA( c ) lexp(/a8 t ) ,  / 

1 2(d , t )  *  SR(d)  +  /87a  SA(d),  and  I 

0  <'  c,d  <  t  <_  tf.  / 


p(a,0,t)  and  p(b,0,t)  are  positive,  monoton ical 1 y  increasing  functions  of  t, 
as  previously  mentioned,  and  thus  SA(t)  and  SR( L )  are  also.  It  is  noted  that 
although  SA(t)  and  Sh(l)  are  functions  of  a  and  8,  they  are  absolutely 
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a0 

bounded  for  any  t  >  0.  This  comes  about  from  the  fact  that  £  p(a,0,«)  *  P(A), 

b0  a"1 

the  probability  A  wins,  and  £  p(0,b,°°)  =  P(B),  the  probability  B  wins. 

b-l 

Therefore,  certainly  S^(t)  <  a  a^PCA)  and  SR(t)  <  B  bQp(B)  where  p(A)  +  P(B)  =  1 
Thus,  by  inspection  of  equations  (12)  and  noting  that  l^Cd.t)  >  0  for  t  >  0, 

(L)  If  either  A  or  A^  is  negative  at  anv  time,  t,  (or  any  range  of 
times)  then  the  other  must  be  positive. 

(2)  If  A  (t)  or  Ag(t)  is  zero  for  some  particular  values  of  a.B.a^.bg,  ami  t, 
then  the  other  mus t  be  positive  at  these  points.  Such  points  are  crossing 
points  for  the  A  with  the  zero  value. 

(3)  Both  differences  may  be  positive  but  from  (L)  above  they  cannot  both  be 
negative  at  the  same  time.  In  particular  they  are  both  positive  for  strict  L 
parity  (i.e.  a^  =  b^,  a  =  8). 

And  from  many  particular  examples  in  Part  Two-l,  we  have 

(4)  A  or  A  may  be  negative,  positive  or  zero  and  Lhese  differences  can  be 
as  high  as  about  30  or  40  percent  of  the  initial  values. 

Furthermo re , 

(5)  Although  all  the  above  applies  for  t  *  l  p ,  in  fact,  it  can  be  sL<wn  they 
are  also  true  for  true  t  >  tp. 

(b)  Finally,  from  (l)  thru  (5)  above,  nontrarv  to  some  statements  in  the 
literature  (for  example,  see  Farrell  (1976,  p.  3)),  the  1.  enuations  cannot  he 
used  as  universal  bounds  on  the  ST,  mean  value  functions. 

2.  The  Linear  Law 


Taylor  (1983,  equation  4.12.24  p.  303)  has  shown  for  the  Linear  Law  that 
it  is  easy  to  derive  the  relation  A^(t)  =  (3/a)A  (t)  which  clearly  indicates 
that,  in  this  case,  the  differences  are  always  of  the  same  sign  and  a 
crossover  in  one  is  accompanied  by  a  crossover  in  the  other  at  exact lv  the 
same  time. 

f.  The  Difference  gear  Tlwe  Zero  for  Certain  Special  Cases 

Many  authors  have  noted  that  for  the  Suuare  Law  (even  with  variations 
such  as  break-points  or  reinforcements)  the  I,  and  ST,  mein  value  differential 
equations  are  identical  if  absorption  probabilities  are  taken  to  he  noligihle, 
that  is  set  them  equal  to  zero.  This  only  makes  sense  for  very  large  initial 
numbers  and  for  early  In  the  combat.  Snow  (1948)  first  mentioned  this  and  it 
was  later  exploited  by  Marshall  (1985),  Clark  (1989)  and  Koopman  (1970). 

Grainger  (1976,  Appendix  G,  pp.  87  and  88)  has  shown  that  for  L  models 
given  by 

x'(t)  »  -6xcv^,  y’(t)  =  -ax'vf,  (13) 
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where  c,d,e,f  are  positive  integers  or  zero,  the  corresponding  SI*  mean  value 
functions  (near  enough  to  t*0  so  that  the  absorption  probabilities  are 
negligible)  are  given  by 


rn^(t)  =  -BE  [ACBd  ] ,  ra^(t)  =  -aE[AeBf  ]. 


(14) 


In  general,  A(t)  and  B(t)  are1 not  independent,  and  therefore  functions  of 
A(t)  and  B(t)  are  not  independent  and  the  only  values  of  c,d,e,  and  f  for 
which  (13)  and  (14)  coincide  (near  t=0,  of  course)  are  (c=0,  d=l),  (c=l,  d= 0) 
and  (c^O.d^O)  for  the  first  equation  along  with  (e*0,  f  =  l),  (e=l,f=0)  and 
(e^O.f^O)  for  the  second  equation.  This  gives  nine  possible  combinations  of 
which  only  (a)  (c^O.d^l ,eal ,f =0),  the  standard  Square  Law  and  (b)  (c=d=e=f=0), 
the  Linear  Law  mode L  where  the  battle  is  a  sequence  of  one-on-one  duels,  are 
interesting  situations.  For  these,  the  difference  is  nearly  zero  in  the 
neighborhood  of  L*0.  For  the  general  Linear  Law  case  (and  other  situations) 
this  indicates  that  there  may  be  large  differences  even  near  time  zero.  Snow 
(1948,  p.  25)  earlier  came  to  the  same  conclusion  for  the  Square  Law. 

G.  Sumary 

Summarizing  the  major  points  in  this  socLion: 

(1)  All  L-SL  mean  value  equivalent  pairs  differ  for  all  times  except  at 
possible  crossing  points.  These  differences  mav  be  very  large. 

(2)  At  least  for  the  Square  Law,  the  I.  trajectories  ire  neither  a  universal 
upper  or  lower  bound  of  the  SL  mean  value  trajectories. 

(3)  F.ven  near  time  zero,  the  L  and  SL  mean  value  trajectories  mav 

differ  considerably  (they  do  not  differ  material lv  for  the  Square  Law  and  the 
sequence  of  one-on-one  duels  version  of  the  Linear  Law). 

III.  FALLACIES  IN  CONTINUOUS  STATE-SPACE  APPROXIMATIONS  FOR  THE  SL  MEAN 

VALUE  FUNCTION 

There  have  been  several  attempts  to  show  that  specific  SL  models  converge 
in  probability  (in  some  sense)  to  L  equivalents.  This  is  another  effort  to 
show  that  L  is  a  good  approximation  to  SL  for  large  numbers.  These  attempts 
have  been  widely  misinterpreted  and  generally  misunderstood.  Tt  should  he 
understood  that  our  analyses  here  reinforce  the  results  in  Section  TT  above. 

A.  Rigorous  Convergence  In  Probability 

Etter  (1971)  and  Karr  (1976)  have  shown  rigorously  that  transformed 
versions  of  the  SL  Linear  and  Square  Laws  converge  in  probability  to  the  L 
laws.  Karr  also  shows  that  transformed  Mixed  (one  side  Linear,  one  side 
Square)  and  Square  Law  versions  with  continuous  reinforcement  also  converge  in 
probability  to  the  corresponding  L  laws.  Their  proofs  differ  (Karr  uses 
probability  arguments  and  Etter  relies  on  function  theory)  but  the 
transformations  are  essentially  the  same.  In  both,  a  state  space  with  Initial 
conditions  (lcs)  of  (aQ.bg)  is  expanded  in  discrete  jumps  (keeping  the  ratio 
of  lcs  constant)  to  lcs  (2aQ,  2b^ ),(  3an,  3Hq  ),...,  (ka^,  kbp,).  In  Karr,  the 

rvs  Ajt(t),  Bjc(t),  where  the  subscript  k  refers  to  the  kL^  space  in  the 
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sequence  above,  are  Chen  transformed  by 

Xk(t)  -  Ak(t)/k,  Yk(t)  =  Bk(t)/k.  (15) 

This  takes  the  expanded  space  with  lattice  points  defined  on  all  non-negative 
integer  pairs  (a,b),  except  (0,0),  such  that  a  <  ka~,  ^  <  kbg  and  transforms 
it  to  a  new  (Xk,Yk)  space  with  ics  (ag,bg)  and  exactly  the  same  number  of 
lattice  points  as  the  untransformed  space  hut  now  they  are  spaced  l/k  units 
apart  in  both  the  x  and  y  directions  (see  Figure  7).  Of  course,  the  mapping 
in  (15)  brings  exactly  the  probabilities  in  the  untransformed  space  onto  the 
corresponding  lattice  points  in  the  transformed  space. 


It  is  very  important  to  observe  that  in  any  transformation  of  the  type 
given  by  equation  (15),  if  a  comparison  is  to  be  made  with  the  corresponding 
Lanchester  equations  the  same  transformations  must  be  made  on  the  L  equations 
(see  Figure  7,  in  which  and  y^  are  the  standard  L  solutions). 

Etter's  transformation  replaces  equation  (15)  by  X^t)  =  A^tK^/k, 

Yk(t)  *  Bk(t)n^/k,  where  Sq.Hq  are  positive  constants  <  l.  He  also  implies 
(correctly)  that  this  transformation  is  similar  to  *k(t)  =  ^(tK^/kta^  +  ) » 

Yk(t)  =  Bk(t)pQ/k(ag  +  b^).  The  corresponding  ics  on  these  two  differently 
transformed  spaces  are  (a0S  ,  b0n0 )  and  (an?r)/(af)  +  b^),  b0n0/(a0  +  b^ ) ) 
respectively  and  the  lattice  celi  sizes  are  £(^/k  bv  n^/k  and  F^/lc(an  +  b^)  by 
n()/k  (a0  +  b^),  respectively. 

There  are  no  essential  differences  in  all  these  transformations  as  Ihev 
all  reauire  that  the  ics  on  the  untransformed  spaces  go  to  infinity  at  the 
same  rate  and  that  the  cells  in  the  transformed  spaces  retain  their  shapes 
(the  ratio  of  the  sides  are  constant)  and  decrease  uniformly.  Thus  as  V  *  <*• 
the  untransformed  space  increases  without  limit  and  the  transformed  space 
remains  the  same  size  and  shape  hut  the  number  of  lattice  points  increases 
with  decreasing  distances  between  them,  (see  Figure  7). 


It  should  he  noted  that  for  some  of  Karr's  results,  he  also  dilates  time 
by  t/k. 


To  clarify  the  situation  let  us  use  the  Karr  transformation  direct l v  on 
the  Square  Law  SL  state  equations.  To  the  best  knowledge  of  the  authors  this 
has  not  been  done  in  this  manner  before.  First  we  define  for  ap.b,-,  and  k  = 

1,2,... 

Pk(a,b,l)  =  PfAk(t)  =■  a,  Bk(t)  =  b,  ic(kaf),khn)  |  , 


=  0,  a  >  ka^  or  b  >  kh^  , 
=  0,  a  <  0  or  b  <  0  , 

■  0,  a  ■  0  and  b  *  0  , 

-  0,  t  <  0  , 


(  lh) 


1,  t  ■  0,  a  -  ka^,  b  *  kb^ 
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SL 


discrete  spaces 


I  -  expanding  the  state  space 
ri  -  transformation  (equation  (!">)) 

[If  -  Limiting  situations  as  k  *  <* 

Figure  7.  Convergence  in  Prohahilitv  for  Fverv  Fixed  Tfme,  t 


Thus , 


3pk(a,h,t) 

- -  =•  'xa(pk(a,b+l  ,t )  -  pk(a,b,t)| 

+  8bfpk(a+l ,b,t)  -  pk(a,b,t)l 

with  ic  p^CkaQ.Icbfl.O)  =  l  (see  Clark  (1969,  p.  69)).  Note  that  for  Lhe 
limiting  process  (where  k  ♦  ®)  only  interior  points  need  be  considered  an(i 
boundary  equations  are  Ignored.  Using  transformation  equations  (15),  which 
imply  that  x  ■  a/k,  y  *  b/k  (where  x  and  y  are  values  of  Xk(t)  and  Y^(l) 
respectively),  and  the  mapping  is  on  Lhe  probabilities  onlv  (i.e.  p^(a, b,t) 
Pk(x,y,t)), 
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(18) 


3pk(x,y,t)  ox[pk(x,y+l/k,t)  -  Pk(x,v,t)] 
3t  "  Uk 

8y[pk(x+l/k,y,t)  -  Pk(x,v,t)] 

+  iTk 


with  ic  pk(a0,bg,0)  =«  l. 

The  probabilities  in  equation  (18)  are  on  all  the  lattice  points  x,y  in 
the  transformed  space,  where  the  lattice  points  are  rational  number  pairs  on 
the  rectangle  { lx [0,b^ |} .  However,  consider  any  fixed  point  (x^,y^), 
rational  or  irrational,  in  the  rectangle.  Each  is  in  some  cell  bounded  bv  the 
lattice  points  for  every  k  (see  Figure  8  below).  The  cell  sizes  are 
diminishing  as  k  Increases.  8y  replacing 


Figure  8.  Transformed  SI.  St  it. ‘-Space  it  Each  I. 

x  by  +  C ,  y  by  y^  +•  n ,  x  +•  1/k  by  +  £  +•  1/k,  and  v  +■  1/k  bv 
y^  +  n  +  l /k  equation  (18)  becomes 

3pk(Xj  +5,yl  -t-n.t)  a(Xj  H  )  (  Pk(Xj +C  ,y  (  +  n  <-  l/k,t)  -  pk(x|+^,v1  +  n  ,  t )  ] 

—  -  — 

8(y1^n)(pk(x1H  +  1/k,  y(+n,t)  -  pk(xj+^,  yt+n,t)l 
+  - m - 

with  ic,  P]t(aQ,bg,  0)  *  1.  Note  that  0  <  n,£  4  1/k,  therefore  as 

k  ♦  «,  n  ,5  ♦  0. 
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Now,  essentially,  Karr  and  Etter  have  shown  that  there  exists  a 
continuous  function,  p,  with  first  derivatives  defined  everywhere  on  the 
transformed  rectangular  state  space  which  p^  tends  to  in  the  limit 
as  k  +  *.  Taking  the  limit  as  k  ♦  *  on  both  3ides  of  (19)  and  dropping  the 
subscripts  on  x  and  y  yields 

3j>(x,,y,t)  „  flx  ££(2!u£a1>  +  gy  .  (20) 

3t  3y  3x 

In  the  limiting  process,  the  p  functions  have  gone  from  a  joint 
probability  mass  function  on  the  Lattice  points  to  a  joint  probability  density 
function  (pdf)  on  the  transformed  continuous  rectangular  state-space.  Thus, 
the  initial  condition  is  given  by, 

p(a0,b0,0)  *  5(x  -  aQ)  S(y  -  bQ).  (21) 

Williams  (1963,  p.  38),  Koopman  (1970,  p.  870)  and  Taylor  (1972,  o.  1-44) 
have  shown  that,  using  the  method  of  characteristics 

dx  dj/  dJL_ 

By  «x  -1  0  ’ 


the  solution  is 


(x,y,t)  =»  constant,  for  x,v  satisfying, \ 
'(l)  =  -By,  y'(t)  =  -ax,  with. 


p( 
x 

x(0) 


a0* 


v(0) 


V 


! 


(121 


At  this  point  it  is  useful  to  examine  the  I.  equations  corresponding  to 
the  expanded  space  SL  equations.  Thev  are 


x'(t)  =  -By,  v'(t)  =  -ax 


with  ics 


x(0 )  »  ka(),  v(  0 )  =»  kb^. 


All  points  on  the  k  3  l  curves  are  transposed  upward  bv  l  factor  of  k. 
Correspond  ing  to  X^(t)  and  Y^(t),  x  and  y  must  he  transformed,  as  follows,  to 
Xj  *  x/k,  y  =  y/k  to  obtain  the  equations  given  above  with  X| ,  Vj  replacing 
x,y  and  with  ics,  Xj(0)  »  aQ  and  y^(0)  =  Sq.  So  the  proper  equations  to 
compare  with  equation  (20)  indeed  are  given  by  equations  (22).  Equations  (22) 
are  the  standard  L  Square  Law  equations  whose  well-known  solution  shall  he 
caLled  Xj(t),  y^(t)  with  x^(0)  =■  aQ  and  y^(0)  =  6q.  This  means  that  at  every 
t  In  the  interval  (0,  tp) 

p(x,y ,  t)  »  5(x-Xj)  (  y  — y  (  > .  (23) 
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Thus,  It  has  been  shown  that  \(t)  converges  in  distribution  to  p(x,t)  - 
6(x-x,),  that  is,  all  the  probability  mass  is  concentrated  at  the  point  x^(t) 
for  every  t  in  (O.t^)  and  similarly 

Yk<t)  X  <S(y-yL>* 


Now,  when  a  rv  converges  in  distribution  to  a  constant  that  Lmplies  the 
stronger  condition  that  it  converges  in  probability  to  the  same  constant,  see 
(Li,  p.  246].  Therefore,  for  0  <  t  <  tf, 

Xk(t)  X  xL(t), 

Vk(t)  X  yf (t ) 


(24) 


Furthermore,  in  the  transformed  space  of  Lhe  rectangle  {  [0,a~  1*  fO,h  ]}  , 
a  statement  concerning  how  expected  values  K[X^(t)]  and  E[Y^(t)]  behave  can  he 
made.  It  comes  from  a  consideration  of  the  convergence  in  probability 
equations  (24)  which,  when  written  in  terras  of  the  definition  of  convergence 
in  probability,  become  (for  the  A  side  only) 

lim  P( |X  (t)  -  x  ( t ) |  <  e]  -  l.  (25) 


The  interpretation  of  equation  (25)  is  that  in  the  limit  as  k  +  °°  the  total 
mass  of  the  sequence  of  random  variables  tX  (t)]  becomes  concentrated  at 
point  x^(t).  Now,  the  X,  (t)  can  have  positive  probability  only  at  discrete 
lattice  points  0,  l/k,  2/k,...,  a0  in  the  hounded  interval  (0,apj.  Thus, 
as  k  ♦  ®,  any  probability  at  these  Lattice  points  (unless  x^( t)  happens  to  he 
a  lattice  point  for  certain  k  values)  approaches  zero.  Therefore,  the  first 
moment  contribution  about  zero  from  any  mass  it  lattice  points  approaches  zero 
(except,  of  course,  if  x, ( t )  happens  to  he  at  i  lattice  point  for  certain 
ks).  Clearly,  then  there  ts  a  verv  important  conclusion  in  the  transformed 
space  that 


lim  E(Xk(t)I  =  x^(t).  (26) 

k 


From  equation  (26)  it  is  clear  that  KfX  (t)|  =  x((t)  *c(k),  where 
lim  e(k)  *  0.  Thus,  taking  expected  values  in  equation  (15)  yields 
k+oo 


or 


F.  ( A  ( t  )  ] 

- - -  -  x^(t)  +  e(k) 


El\(t)l 

kxL(t) 


l  + 


e(k) 

xL(t) 
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and 


IfA^Ct)] 


lira  ,  f  r\ 
lc~  k  xL(t) 


E[B  (t)l 

lira  - - TTv 

k>«  k  yL(t) 


Similarly  for  Y^(t), 


(27) 


Equations  (24),  (26)  and  (27)  are  the  principal  focus  of  this  section. 
(24)  and  (26)  are  a  Law-o f -Large-Numbers  type  result  for  the  Lanchester  Square 
Law.  It  is  important  to  note  that  these  two  equations  say  that  the 
transformed  SL  rvs  converge  in  probability  to  the  L  Square  Law  which  is  the 
limiting  mean.  Contrary  to  many  statements  in  the  literature  (both  explicit 
and  implicit)  this  does  not  say  there  is  convergence  of  the  means  (i.e.  mean 
SL  may  not  tend  to  L)  in  the  untransformed  space.  All  that  can  be  said  about 
the  untrana formed  SL  sequence  is  given  in  equations  (27).  It  has  been  noted 
previously  that  for  A^(t),  \(t)  that  the  corresponding  L  results  are  kx^(t) 
and  ky^(t).  Thus,  equation  (27)  says  that  the  ratio  of  expected  values  of  SL 
to  corresponding  L  equations  go  to  a  limiting  value  of  l.  This  does  not 
necessarily  mean  that  [E(A-  (t)J  -  kXj(t)J  goes  to  zero  in  the  limit.  .As  a 
matter  of  fact,  this  difference  may  go  to  a  constant  (Including  zero)  or 
infinity  and  the  ratio  still  go  to  l.  The  authors  believe  that  this  important 
distinction  is  pointed  out  here  for  the  first  time  in  Lanchester  literature. 
Mothing  that  has  been  done  to  date  says  anything  more  about  the  limiting 
behavior  of  the  untransformed  difference. 


Equation  (20)  has  been  derived  in  several  other  less  rigorous,  intuitive 
ways  which  seem  to  have  concealed  its  true  message  as  given  above.  These 
shall  now  be  briefly  examined. 

B.  Diffusion  Approximations 

Equation  (20)  can  be  considered  a  first  order  diffusion  theorv 
approximation  to  the  SL  Square  Law  process  and  it  (or  its  implications)  have 
been  arrived  at  in  several  nonrigorous  wavs. 

1.  Taylor 'a  Series  Expansions 

First  is  the  Taylor's  series  approximation  (see  [7|).  Again  the  SL 
Square  Law  is  used  to  illustrate.  The  notion  involved  here  is  to  replace  the 
discrete  function  given  by  equation  (17)  (with  k=l),  by  a  continuous  function 
that  goes  approximately  through  each  of  the  discrete  values  of  p(a,b,t)  on  the 
a,b  axes  at  each  time  t.  Thus,  replacing  discrete  a,b  by  continuous  x,v 
yields 

3  ^*3  t  *  ox[p(x,y+I,t)  -  p(  x  ,y ,  t )  | 

+  8y[p(x+l ,y,t)  -  p(x,y,t)I  . 

Now  expand  p(x,y+l,t)  in  a  Taylor's  series  in  powers  of  l  around  y  and 
p(x+l,y,t)  in  powers  of  1  around  x  and  retain  only  the  first  order  terms  to 
get  equation  (20)  immediately. 
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The  usual  explanation  of  when  this  (l.e.  equation  (20))  is  a  fair 
approximation  to  equation  (17)  is  for  axt  and  Byt  to  be  large.  This  is  only 
speculation  based  on  particular  calculations.  This  really  means  for 
fixed  <*,8,t  that  x  and  y  must  -*■  ®,  which  implies  a  transformation  of  the 
type  previously  discussed.  Therefore,  the  usual  assumption  that  this  applies 
to  the  untransformed  space  is  not  correct.  As  was  shown  earlier  the  only 
correct  conclusion  on  the  untransformed  space  Is  given  by  equation  (27). 

2.  Approximating  Differences  with  Derivatives 

Earlier  Willard  (1962  pp.  31-33)  aLso  arrived  at  equation  (20)  by  a 
somewhat  more  intuitive  approach  where  he  divided  the  first  term  on  the  rhs  of 
equation  (17)  by  (b+l)  -  b  and  the  second  by  (a+1)  -  a  and  called  these 
Ab  and  Aa,  respectively.  He  then  replaced  b+l  by  b  +  Ab  and  a  +  I  by 
a  +  Aa  and  then  assumed  a  and  b  were  continuous  and  arrived  at  equation  (20) 
by  letting  Aa,  Ab  +  0.  This,  of  course,  is  reasonable  only  for  large  a,b  and 
again  amounts  to  a  mapping  because  Aa  and  Ab  are  always  exactly  1  and  can  only 
be  made  a  truly  variable  difference  by  a  transformation.  His  solution  is 
arrived  at  In  a  different  manner  but  is  (as  It  should  be)  equation  (23). 
However,  he  incorrectly  states  that  this  implies  that  L  is  the  limiting 
solution  to  untransformed  SL.  Again  we  reiterate  the  only  conclusion  on  the 
basic  SL  equation  is  given  by  equation  (27). 

Koopraan  (  1970,  p.  8’’  ■  ind  Taylor  (1972,  p.  1-42)  also  used  this 
technique  to  get  equation  ,20)  and  solved  it  by  the  method  of  characteristics. 

Helmbold  (1966.  pp.  632-635)  also  uses  this  technique  for  discrete  state- 
space  and  discrete  time  parameter  models  of  marksmen  versus  passive  targets 
and  many  versus  many  battles  to  get  L  equations.  This  involves  simultaneously 
passing  to  derivatives  from  differences  on  both  state-space  variables  .and 
time.  The  discrete  time  parameter  in  these  models  comes  about  because  all 
contestants  on  a  side  fire  in  volleys  .at  discrete  time  Intervals. 

3.  Tl«e-Independent  State-Space  Analyses 

Williams  (1963,  p.  31,  et.  seq. )  has  written  difference  equations  on  the 
moments  of  the  terminal  survivor  distribution  as  a  function  of  the  initial 
conditions  for  both  the  Square  Law  and  the  Linear  Law.  This  is  done  bv  a 
random  walk  on  the  time— independent  joint  state-space.  For  example,  he  shows 
that  for  the  Square  Law 


8  b. 


VvV 


oa0+6ba 


aa0 

Uk(an“l’b0)  +  aaf)+6  b 


\(Vbo_,) 


(29) 


where,  U,  is  the  k^—  moment  about  the  origin  of  the  marginal  distribution  of 
the  A  side  survivors  (l.e.,  when  A  wins).  He  uses  the  Taylor's  series 
expansion  In  powers  of  one  (as  explicated  earlier  in  a  different  context)  and 
again  retains  only  first  order  terms  to  get  the  diffusion  expression 


8  b. 


3uk(a0,b0) 

3a„ 


+  aa. 


3uk(a0,b0) 


3  b, 


=  0. 
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It  Is  recognized,  of  course,  that  this  really  represents  an  equation  for  a 
Karr-type  transformation  and  that  ag  and  bg  should  he  replaced  by  variables  on 
the  transformed  space.  He  then  solves  this  to  show  that  urt  ■  Prob  [A  wins]  =  1, 

(~~2  2  2  U 

for  all  real  u^a^bg)  =  /  a^  -  (S/a)  b^  and  n,  =  U.;  therefore  the  variance 
is  zero  and  furthermore  all  Higher  moments  around  the  mean  are  zero.  These 
equations  are  valid  only  for  ag/ a  >  bg/0  of  course.  Similar  expressions  can 
be  arrived  at  for  the  case  where  B  wins.  This  merely  shows  that  on  the 
transformed  space  the  terminal  distribution  is  the  terminal  Lanchester 
point  as  one  expects  from  the  Karr-Etter  development. 

He  shows  the  same  thing  for  the  Linear  Law. 

Covey  (1969  p.  31)  has  shown  that  for  the  Linear  Law,  Am^(t)/A'<  =  -S/(a-H3), 
where  k.  is  the  number  of  the  transition  on  the  macgtnal  time-independent 
sLate-space  points.  This  is  exactly  true  and  shows  that  the  marginal  mean 
decreases  linearly  on  equally  spaced  marginal  points  that  lie  exactly  on  the  L 
Linear  Law  solution  in  state-space.  However,  at  this  point  he  makes  the 
grossly  incorrect  inference  that  (for  large  numbers)  the  expression  above  may 
he  approximated  by  dmA(t)/dt  3  -6/(a+g)  which  is,  of  course,  not  true  since 
it  implies  jumping  from  time-independent  state-space  to  the  time  trace  which 
is  totally  unjustified  and,  of  course,  does  not  give  Lhe  correct  result. 

4.  Miscellaneous  Analyses 

(1)  Gye  and  Lewis  (1974  pp.  6-7)  give  a  curious  twist  to  all  this  Sv 
starting  with  Lhe  usual  spurious  notion  that  equation  (20)  applies  to  the 
untransformed  space  and  then  applying  the  fundamental  calculus  identity, 

dp/dt  =  3p/dt  +  (3  p/3 x)(dx/dt)  +•  ( 3 p/3 y  )(dy/d t)  to  arrive  at  the  conclusion  that 
the  Square  Law  Lanchester  equations  are  really  Lhe  mo da  1  trace  rather  than  the 
expected  value  trace  in  the  untrans formed  space.  It  has  been  shown  earlier  that 
it  is  not  the  expected  value  trace  and  this  certainly  does  not  establish  it  as 
the  modal  trace  as  equation  (20)  only  applies  to  the  transformed  space  in  Lhe 
limit  and  no  meaning  can  be  attached  to  Lhe  expressions  dx/dt  and  dv/dt  in  the 
untransformed  space.  In  the  transformed  space  it  has  been  shown  that  all 
probability  is  located  on  Lhe  Lanchester  trace  and,  therefore,  all  moments, 
etc.  are  located  there  and  to  say  that  the  mode  is  located  there  adds  nothing. 

(2)  Farrell  (1976,  pp.  4-13)  has  made  an  attempt  at  bounding  the  solutions 
to  the  mean  value  functions  for  the  SL  Square  Law.  Tn  Lhe  first  part  of  his 
paper  he  erroneously  states  that  the  L  equations  are  lower  bounds  on  the  SL  mean 
value  functions  (mA(t)  and  mB(t))  for  all  an,b0,  t  <  Lf,  a, 6.  This  is  patently 
not  true  from  the  discussion  in  Section  TT  and  from  numerous  counter  examples. 

He  then  extends  this  to  the  Linear  Law  by  some  assumptions  which  are  prohahlv 
true  but  leave  him  with  the  same  defect  that  vitiated  his  conclusions  in  Lhe 
Square  Law.  Farrell's  argument  is  based  on  the  fact  that  equations  (2)  for  the 
SL  mean  values  would  look  exactly  like  the  Lanchester  equations  if  It  were  not 
for  the  two  nonnegative  terms  on  the  right  hand  sides  involving  the  absorption 
probabilities'.  For  any  positive  t,  these  prohabl  1  i  t  les  become  positive  and  one 
might  quickly  conclude  that  the  SL  solutions  are  upper  bounds  to  the  L  survivor 
functions.  However,  because  of  the  interaction  between  the  equations  of  each 
pair  this  conclusion  is  not  true;  and  in  fact,  there  are  manv  counter  examples. 
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In  the  next  part  he  attempts  to  develop  upper  bounds  on  these  same 
functions.  Again  there  is  a  flaw  which  vitiates  the  results.  He  makes  a 
transformation  which  does  not  require  the  initial  force  sizes  to  go 
to  infinity  at  the  same  rate  and  thereby  destroys  the  validity  of  the  work. 
This  can  be  demonstrated  as  follows.  The  Farrell  transformation  is  given  as 

X(t)  *  A(t)/a^  and  Y(t)  =  B(t)/b^. 

When  this  transformation  is  applied  to  equation  (17)  with  k  =  l,  the  result- is 


5 1 =  «xaQ[p(x,y  +  l/b^.t)  -  p(x,y,t)] 


+  8ybQ[p(x  +  I/a^.v.t)  -  p(x,y,t). 


ax 


aQ  [_p(  x,y  +  l/b0>t)  -  p(x,v,l) 


"0 


1  /  h 


0 


b„  [p(x  +  l /sq.v, t)  -  p(x,y,t)' 

+  8v - P7 - 

a0  lM0 


How  let  ag/b0  =  a  constant  say,  c,  and  let  an  and  b^  *  "  (using  the  same 
arguments  about  cells  on  the  transformed  state  space  as  before)  the  result  is 


1  p  ( x ,  v ,  t ) 

3l 


axe 


3p(x,y,t)  + 

•1  V 


i  3p(x,v,t) 

C  3  X 


no) 


Now,  as  long  as  c  is  held  fixed  this  transformation 
results  as  before  because  the  corresponding  transformed 


will  give  the  same 
1.  equations  are 


with  ics 


dv/dt  =  -acx,  dx/dl  =  -(5/c)v 


x(0)  =  v( 0 )  =  1, 


which  indeed  give  the  correct  x^  ,  coordinates  for  the  solution  to  equation  (30) 
which  is  again  p(x,y,t)  =  <5  ( x  -  x.)  6(v  -  ).  However,  in  all  his  conclusions, 

Farrell  lets  one  side  have  a  fixed  Initial  condition  and  lets  the  other  sides's  ic 
-*■  -*>.  This  gives  either  c  =  0  or  c  -  0D  and  the  whole  analysis  collapses. 


In  fact,  if  we  let  a^  00  then  c  *  •»  and  the  transition  probabilities 
downward  will  tend  to  one  and  leftward  to  zero  and  in  the  limit  all  probability 
mass  will  be  concentrated  at  p(co,0,t),  for  all  l  for  boLh  spaces. 

The  basic  point  is  that  in  any  transformation  of  this  tvpe,  if  the  initial 
numbers  on  one  side  go  to  infinity  both  must  do  so  and  furthermore  both  must 
do  so  at  the  same  rate. 


IV.  FALLACIES  IN  THE  APPROXIMATION  Of  GENERAL  RENEWAL  PROCESSES  BY 

SL  PROCESSES 


In  all  that  has  been  said  up  to  this  point  it  has  been  assumed  in  the  SL  model 
that  all  interfiring  times  XA,  Xg  are  negative  exponent i a  I Iv  distributed  (ned) 
rvs.  This,  of  course,  greatly  simplifies  matters  as  the  ned  rv  has  "no  memory". 
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Two  attempts  have  been  made  to  justify  using  the  ned  SL  model  as  a  good 
approximation  to  the  general  renewal  model  (l.e.  with  general  ifts,  X^,  Xg),  under 
certain  circumstances.  These  two  arguments  are  now  examined  in  some  detail. 

A.  The  Individual  Firer  Argument 

First,  there  is  the  individual  firer  argument  which  is  widely  used  in  this 
country  and  is  based  (see  Bonder  and  Farrell  (1970),  pp.  84-86)  on  Blackwell's 
Theorem  in  renewal  theory.  This  argument  goes  as  follows: 

(1)  Observe  each  independent  firer  on  one  side,  all  with  identical 
independent  ifts  which  are  general  rvs ,  X  with  mean  Q  and  with  kill 
probability  p.  Embedded  in  each  of  these  processes  is  a  killing  orocess  with 
mean  C/p,  which  is  also  iid  and  therefore  renewal  theory  is  applicable. 

Set  s /p  -  u ,  for  simplicity,  and  call  the  embedded  stochastic  process  Z(t). 

(2)  Each  firer  has  a  sequence  of  kitLs  which  are  counted  and  which  is  a  rv 
in  time,  X(t),  the  counting  distribution  of  Z(l). 

(3)  The  uncountably  infinite  ensemble  of  all  realizations  of  M( t )  has  a  mean 
value  function  of  time,  n(t),  whose  slope  -*■  1/u  as  t  -*■ 

(4)  (3)  implies  that  Z(t)  *  ned  with  mean  u,  which  is  not  true  (as  shall  he 
shown ) . 

(5)  (4)  implies  that  a  superposition  of  all  firers  on  one  side  tends  to  have 
ned  interkilling  times.  if  (4)  is  not  true,  this  is,  of  course,  not  true. 

Figure  9  below  shows  a  few  of  the  uncountablv  infinite  set  of  possible 

realizations  of  N(t).  Of  course,  at  overv  time  L , _ L  he  mean  is  the  average 

over  the  entire  ensemble  and  is  denoted  K(N(t)|  =  nit). 

In  Figure  10  we  graphically  show  some  definitions  from  ordinarv  (non- 
terminating)  renewal  theory  and  below  are  some  theorems  from  ordinarv  renewal 
theory  (see  reference  f  8  )  ) : 

(1)  Lim  n(t)/t  =  l/u  (ElemenLarv  Renewal  Theorem),  this  means  the  time 
average  (slope  of  the  chord)  of  the  mean  value  function  tends  to  I /u . 

(2)  Lim  N(t)/t  =  l/u  with  probability  one,  this  means  that  the  time  average 
(chord)  of  every  evolution  tends  to  l/u  (not  shown  in  Figure  10,  see  Figure  9). 

(3)  Lim  dn(t)/dt  =  l/u  (Blackwell's  The<  >rem  or  the  Key  Renewal  Theorem) 

t+“ 

means  that  the  instantaneous  rate  tends  to  l/u. 

It  is  (3)  above  which  has  been  invoked  to  justify  using  P./u,  for  n  and 
Pg/Ug  for  8  in  SL  and  L  approximations  to  OR. 

The  fallacies  in  this  assumption  are  several; 
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(2)  n(t)  is  an  ensemble  average  and  therefore  rates  defined  in  properties  (I) 
and  (3)  above  are  ensemble  rates.  Thev  are  a  weighted  average  of  all  rates  for 
certain  special  mutually  exclusive  and  exhaustive  subsets;  weighted,  as  we  shall 
see,  by  the  probability  that  the  backward  recurrence  time  is  some  particular 
value.  Thus,  they  cannot  he  used  indiscriminately  in  probability  calculations, 

(3)  none  of  these  theorems  apply  to  terminating  processes  (which  are  dealt  with 

here),  and  thus  none  of  this  implies  that  the  embedded  killing  process  is 
tending  to  ned  with  mean  In  fact,  for  the  classical  terminating  process  it 

is  easy  to  show  that  lim  dn/dt  =  0. 

To  further  illuminate  point  (2)  above,  examine  what  the  correct 
instantaneous  rate  is.  The  first  thing  to  note  is  that  except  for  ned 
tnterevent  (interkilling)  times  it  is  not  sufficient  to  simply  specify  n(tl 
for  the  state.  To  completely  and  unambigously  define  the  state,  n(t)  and  the 
backward  recurrence  time  (y)  must  be  specified  (see  Figure  11).  This  is 
because,  in  general,  the  system  has  Markovian  "memory"  and  it  does  "remember" 
the  time  of  the  last  event. 


Figure  11.  The  Backward  Recurrence  Time 

The  rv  Y  with  value  y  at  t,  as  shown,  is  the  necessary  additional 
information  to  completely  specify  the  state.  A  function  of  v,  riv)  is  the 
ensemble  rate  for  the  subset  of  all  realizations  with  Y  =  v  and  can  he  used 
for  probability  calculations,  l.e.  r(y)A  =  '’[event  (kill)  in  ( t ,  t  t-A  ) )  where 
r(y)  A  f^(v)/F^(y)  which  is  the  Instantaneous  rite  given  Y  =  v  and  is  the 
weighting  factor  referred  to  In  (2)  above..  Note  that  if  Z  is  ned  Lhen 
r(y)  =*  1 /u  and  does  not  depend  on  y. 

B.  The  Superposition  Argument 

Now  examine  the  second  argument  which  is  called  the  superposition 
argument  and  which  is  widely  used  In  the  Russian  literature,  VentseL  (l°f>4). 
This  argument  bolls  downs  to  the  Pa lm-Khintchine  Lheorem,  see  reference  [9), 
which  essentially  states  that  if  all  Hq  lid  tnterkilling  processes,  Z^,  on  sav 
the  A  side,  are  superposed  to  form  a  new  Interkilling  process,  W^(l), 
then  W^(t)  ♦  ned  with  mean  l/c,  as  a^  ♦  «  if 

ao 

( l )  l  l  /u A  *  c,  and 
i-1  * 
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(2 )  each  uA  -*■  ®, 

where  the  uAs  are  the  means  of  the  ltd  ZAs. 

In  words  this  simply  means  that  for  either  side  (say  the  A  side)  the 
superposed  embedded  killing  process  tends  to  ned  as  the  number  of  simultaneous 
firers  tend  to  infinity  if  each  firer’s  inLerkill  time  mean  tends  to 
infinity.  No  one  has  investigated  how  large  must  be  and . how  large  u ^  must 
be  for  this  to  be  practically  useful. 

It  should  also  be  noted  that  this  theorem  applies  to  non-terminating 
processes  (no  theorem  like  this  for  terminating  processes  is  known  to  the 
authors).  Indeed,  obviously,  as  Lime  progresses  in  a  terminating  situation 
the  requirement  for  large  numbers  will  sooner  or  later  be  badly  violated. 

It  should  be  noted  here  that  the  only  free  world  example  of  this 
misapplication  of  the  Palm-Khintchine  Theorem  that  we  have  found  is  in  Cho 
(1984)  where  he  incorrectly  assumes  it  Is  applicable  in  a  multiple  marksmen 
versus  multiple  passive  targets  situation. 

C.  Firing  and  Killing  Rates 

Next,  an  interesting  controversy  which  continues  to  crop  up  from  time  to 
time  (see  reference  [4])  is  considered.  It  is  the  question  of  Lhe  appropriate 
measure  for  the  individual  firing  rate  (this  is  equivalent  to  the  question  of 
the  measure  on  the  individual  kill  rale  as  it  is  easy  to  prove  that  in  the  f!R 
model,  say  for  the  A  side,  the  individual  kill,  rate  is  p^rA>  if  rA  =  l  /u  is 
defined  as  the  individual  firing  rate).  The  question  which  has  been  raised 
is,  what  Is  correct,  r^  =  1 / K  f  X  ^  ]  or  r^  =  l\fl/X^|7  bonder  [5|  originally 
proposed  that  Rfl/X.]  was  the  proper  rate  hut  Bar foot  [2]  made  an  intuitive 
argument  that  L /K  f  X « ]  Is  the  correct  one  and  implicitly  (though  not 
explicitly)  Bonder  |6]  finally  agreed. 

The  fact  is  that,  in  general,  as  has  been  shown  above,  neither  is  correct 
and  there  is  no  such  general  fixed  rate.  however,  if  one  must  use  such  an 
approximation  there  is  no  question  that  1/F.jX^]  is  the  better  one.  Renewal 
Theory  (see  above)  shows  that  for  Lhe  SI.  process  it  is  exact  and  is  therefore 
a  counter-example  to  the  original  Bonder  thesis;  for  Lhe  (1R  process  it  is 
asymptotically  correct.  From  an  intuitive  viewpoint  what  Lhe  controversy 
bolls  down  to  is  the  following;  suppose  one  collects  n  independent  inlerevoul 
times  ,  l  =  l,2,...n.  Should  one  consider  the  sample  event  rate  to  he 

n/.E.x,  =  l/f.Z.  x,/nl  =  1/KfX.J  or  should  one  consider  each  1/x,  to  He  i 
l“l  l  i=l  l  '  A  n  i 

sample  of  the  event  rate  leading  to  the  sample  event  rate  (  l/x^)/n  -  M ( l / X ^ 

Barfoot  pointed  out  that  the  first  is  an  event  average  and  the  second  is  a 
time  average  and  the  first  is  what  we  seek.  We  note,  in  passing,  that  for  a 
positive  rv,  X,  there  Is  a  theorem  that  F[l/X|  >  1/  K f X |  (see  reference  fill, 
p.  16b).  It  Is  time  to  bury  this  controversy. 

D.  Farther  Coanents 


(1)  It  Is  Important  to  note  that  there  are  iheorems  which  state  that,  in 
general,  even  for  non-termlnatlng  processes,  Lhe  inierevent  times  for 
superposed  renewal  processes  are 
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(a)  not  identically  distributed  and, 

(b)  are  not  independent, 

both  of  ’which  are  necessary  for  renewal  theory  to  apply  to  any  process.  The 
only  exception  is  for  ned  processes. 

(2)  It  should  not  be  implied  that  if,  in  either  the  L  or  SL  formulation,  the 

kill  rates  a,  6  are  made  functions  of  time  that  the  resulting  situation  may 
drive  these  formulations  closer  to  the  OR  model.  What  happens  if  this  is  done 
is  that  in  the  SL  formulation  each  individual  firer's  prodess  becomes  a  non- 
homogeneous  Poisson  process  which  is  indeed  Markovian  but  renewal  theory  does 
not  apply  as  the  interkill  rv  depends  on  time  (i.e.  they  are  not  identically 

distributed)  and  successive  interkill  times  are  dependent.  Thus  the  renewal 

character  of  the  individual  processes  of  the  model  are  destroyed  and  although 

it  is  a  model  of  some  process  it  is  not  necessarily  driving  SL  closer  to  (SR. 

The^  only  possible  function  which  might  have  some  merit  in  this  regard  is 
using  dnA(t)/dt  for  a  and  dng(t)/dt  for  fi,  for  very  large  3q,  b0.  This  has 
not  been  investigated.  The  notion  here  is  that  for  sufficiently  large  a^ ,  b^, 
in  the  superposed  process  there  may  he  enough  backward  recurrence  times  of 
various  sizes  and  in  appropriate  proportions  to  Imitate  the  ensemble  mentioned 
earlier  at  every  t.  Again  though,  the  imitation  will  surely  become  poor  as 
tLme  Increases. 

(3)  A  tvpical  result  from  C.afarian  and  Ancker  (1084)  is  given  fn  Figure  12 
which  shows  a  wide  discrepancy  In  comparable  L,  ST,,  and  OR  model  mean  value 
functions.  This  is  a  two-on-one  situation  with  Frlang  (2)  ifls  on  the  A  side 
and  ned  on  the  B  (one)  side. 

V.  THE  IMPORTANCE  Of  PROCESS  VARIANCE 

Up  untiL  now,  the  focus  has  been  upon  the  notion  that  certain  combat 
processes  are  adequately  described  bv  their  mean  value  functions  and  in 
particular,  hv  a  determini  si ic  Lanchester  approximation  to  this  function. 

A.  The  Initial  SL  Variance 

Implicit  in  the  idea  mentioned  above  is  the  assumption  that  fat  least  for 
large  numbers  on  each  side)  the  variance  of  the  process  is  unimportant.  This 
assumption  is  largely  based  on  Brooks  (10f>3).  Brooks  has  looked  at  processes 
such  as  ours  In  state-space  and  examined  them  at  the  lattice  points  only  and 
in  particular  at  the  lattice  points  where  successivelv  t  event  onlv  has 
occurred,  2  events  only,  3  events  onlv,  etc.  up  to  the  point  where  the  number 
of  events  (k)  is  equal  to  or  less  than  the  min  ( nf) ,  h^ )  (see  Figure  4.)  This 
ensures  that  not  enough  events  (events  being  kills  on  either  side)  have 
occurred  for  an  absorption  to  have  occurred.  This  essentially  means  verv 
early  in  the  combat.  He  then  defines  a  concept  called  "stochastic 
determinism"  as  the  property  of  a  process  at  the  kL^  event,  that  the  quantity. 
Is  small  for  all  k  <  min  and  where  o,  is  the  standard 

deviation  of  the  losses  on  the  A  sine  (this  is  also  the  standard  deviation  of 
the  survivors). 

Next  he  shows  for  the  SL  Linear  [.aw  (in  general,  even  though  lie  claims  it 
for  a  succession  of  one-on-one  duels  only)  that  o^/a^  <  1/2  / a^  ,  k  <  min 
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Figure  12.  An  Example  of  Comparible  h,  SI,,  anil  ;;!<  Model  Mean  Value  Functions. 

(sq,  b(-j),  which  Is  sufficent  for  his  he  f  i  ni  t  ion  of  stochastic  determinism. 
However,  examination  of  the  coefficient  of  variation  of  the  losses  for  Lhis 
case  yields  /a  /kB  .  or  a  minimum  valne  of  /1/mi  n  fa  (. ,  h^ )  /a  IS  .  Clearlv, 
the  term  /a/8  can  be  verv  large  indeed  and  one  would  surely  not  he  justified 
in  calling  Lhis  quantity  small  for  al I  poss ihl e  parameter  values.  Or,  just 
look  at_the  value  of  o^,  which  is  /ka8/(a+8)  tn<l  which  has  a  maximum  of 
(1/2)  /k  .  Attain  for  Large  values  of  a;-y  and  h^ ,  Lhis  can  attain  a  1  arge 
value  since  k  can  equal  min  (af),h(j). 

For  the  Square  haw  Brook's  results  are  similar  (he  only  considers  the 
case  where  a  =8). 

Clark  (1969,  pp.  112-111)  using  a  technique  suggested  lav  Snow  (1949)  has 
given  the  variance  for  the  Sh  Square  haw.  The  technique  solves  the  variance 
equations  again  using  the  assumption  that  the  absorption  probabilities  are 
zero.  This  is  essentially  the  same  idea  as  Brooks  above,  hut  gives  Lhe  more 
informative  time  trace  (good  only  for  times  near  wro).  The  equations  clearlv 
indicate  that  although  the  percentage  losses,  for  large  a^  and  may  he 
small,  the  absolute  values  may  be  large.  This  supports  Lhe  analysis  above. 

Willis  (1982,  p.  6)  has  arrived  at  exactly  the  same  ordinary  differential 
equations  on  the  moments  as  Clark  (1969)  above.  His  technique  is  to  use  the 
differential-difference  equation  (17)  for  interior  points  to  obtain  a  partial 
differential  equation  on  the  joint  moment  generating  function  of  the 
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process.  The  moment  generating  function  is  then  expanded  in  a  Taylor's  series 
in  powers  of  the  two  transform  variables  and  coefficients  of  like  powers  on 
each  side  of  the  equation  are  equated  to  generate  the  ordinary  differential 
equations.  Since  boundary  equations  are  ignored  this  again  implies  absorption 
probabilities  are  zero  and  is  good  only  near  time  zero.  For  some  models  this 
technique  may  be  easier  than  Clark’s  and  Snow’s. 

Weale  (1972  pp.  11-12)  reproduces  Clark's  results  (apparently 
unknowingly).  However,  he  does  give  an  approximate  expression  based  on  a 
Normal  approximation  to* the  joint  pmf  of  A(t)  and  B(t)  for  the  time  (measured 
from  t  =  0)  during  which  Clark's  expressions  are  good  to  any  desired  degree  of 
approximation.  (The  relevant  equations  are  (42),  (43)  p.  20,  (38)  and  (39)  p. 
18  and,  (19)  p.  8  in  Weale.)  The  generaL  idea  is  illustrated  in  Fig  13 
below.  Although  Weale  does  this  for  a  breakpoint  analysis,  the  annihilation 
situation  is  simpler  and  illustrates  the  point.  If  1  -  exp(-cv2)  is  the 
joint  survival  probability  (for  a^,  b~  -*■  <”)  contained  inside  the  contour  c  at 
some  particular  time  for  some  particular  probability  then  certainly  some  of 
the  remaining  probability  is  in  the  form  of  absorption  probability  on  the  a 
and  b  axes.  Roughly  if  a  total  absorption  probability  is  selected  that  is  not 
so  "large"  as  to  distort  the  results  then,  one  minus  this  probability  is  the 
probability  desired  inside  the  contour.  Then,  by  a  trial  and  error  procedure 
the  time  at  which  the  desired  contour  is  tangent  to 


Figure  13.  Survival  Probability  Contours 

one  or  both  of  the  axes  Is  determined.  In  any  event  the  time  Interval  thus 
calculated  is  conservative. 

Perla  and  Lehoczky  (1977  pp.  5-12)  derive  a  diffusion  model  for  the  SL 
Square  Law.  In  this  approximation  the  mean  value  is  assumed  to  be  the  Xj(t), 
yL(t)  time-trace  and  the  pdfs  at  every  time,  t,  are  continous  and  assumed 
Normal  (invoking  the  Central  Limit  Theorem).  Then  the  procedure  derives  the 
second  moments.  The  variances  are  exactly  as  given  by  Clark  (1969)  above. 
There  is  also  a  covariance  determined  which  Clark  could  have,  but,  did  not, 
give.  The  Perla  and  Lehoczky  model  is  only  good  until  the  absorption 
probabilities  are  significant,  and  is  essentially  no  different  than  Clark 
(1969).  Again  they  show  that  early  variances  may  be  large. 

In  general,  the  use  of  a  continuous  state-space  implies  a  limiting 
mapping  (good  only  on  a  tranforraed  space)  as  explained  earlier.  But,  for 
approximations,  these  expressions  with  very  large  initial  numbers  are  used  in 
the  untransformed  state-space.  This  is  justified  by  the  following  reasoning 
(see  Figure  14  below). 
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Figure  14.  The  Diffusion  Approximation 

The  figure  only  shows  a  univariate  situation  hut  the  reasoning  carries 
over  to  the  bivariate  situation.  The  probability  mass  (such  as  at  points  n, 
n+1,  n+2)  is  converted  to  a  rectangular  probability  density  around  each  mass 
point,  as  shown,  with  the  same  height,  since  Lhe  intervals  are  exactly  one 
unit  apart.  The  histogram  thus  created  is  approximated  bv  a  curve  which 
nearly  goes  through  the  mass  points  as  shown  (it  only  goes  through  them 
exactlv  in  the  limiting  transformed  space).  The  mean  value  function  is 
approximately  the  L  equation  with  possibly  large  variance  for  very  large 
initial  condition  and  early  in  the  process. 

The  conclusion  is  that  "stochastic  determinism''  as  defined  by  Brooks  does 
not  ensure  that  there  are  not  large  absolute  variations,  even  in  the  early 
stages  of  SL  models  with  large  initial  numbers  on  each  side. 

B.  The  Terminal  SL  Variance 

It  can  be  deduced  from  Weiss  ( 1 9b )  )  that  for  Lhe  SI,  Linear  Law  the 
terminal  distribution  is  asymptotically  normal.  Also  five  and  Lewis  (1 97b) 
show  asymptotic  normality  for  the  terminal  distribution  of  Lhe  SL  Smiare  t,aw. 

Gye  and  Lewis  (  1974,  p.  19)  have  shown  that  for  Lhe  Snuare  Law  Lhe 
terminal  distribution  of  A  side  survivors  for  large  initial  numbers  is 
approximately  Normal  with  a  standard  deviation  of  about  .51  /a  if  the  A  side 

M' ,  ° 

has  overwhelming  superiority  and  about  .  17  a()  for  an  =  b().  In  both 
cases  a  =  8  =  I. 

Watson  (197b)  uses  Martingale  Theory  to  arrive  at  similar  terminal 
results.  This  Involves  a  transformation  of  A(t)  and  B( t)  to  get  a  Martingale 
rv  whose  terminal  properties  are  easily  arrived  at.  However,  the  inversion  to 
get  terminal  properties  for  A(t)  and  B(t)  is  far  from  simple  and  it  Is  not 
obvious  that  there  is  any  computational  gain. 

It  Is  clear  from  many  examples  that  terminal  variances  for  initial 
conditions  of  any  size  are  always  substantial. 
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C.  The  Transient  SL  Variance 


Taylor  (1972),  (1-45,46)  develops  a  second  order  diffusion  approximation 
to  the  SL  Square  Law  by  standard  diffusion  techniques.  This  partial 
differential  equation  has  unknown  coefficients;  however,  the  following 
technique  identifies  the  coef ficients. 

The  procedure  is  the  same  as  given  previously  starting  with  equation  (28) 
and  expanding  the  p(x,y,t)s  in  powers  of  one  around  x  and  v,  except  that  this 
time  one  keeps  second  order  terras.  It  is  simple  to  show  that  one  immediately 
gets 

=  ax  3_Pi Mdl  +  v  l P<x2vAt) 

3 1  3  v  3  x 

+  ».  (ax  +  By  iJei Spill)  , 

3y  "  3  x 

with  initial  condition 

p(x,v,0)  =  (x-xfl)  ( v — Vq ) . 

This  equation  has  not  been  solved  at  this  time,  and  so  adds  little  to  our 
knowledge.  One  is  tempted  to  hope  that  its  solution  mighL  he  a  Normal  pdf 
with  mean  Xj(t),  v^(t)  and  variances  and  covariance  as  given  by  Clark  (1969) 
and  Per  la  and  Lehockzy  (  1977  ),  however,  ibis  does  not  appear  to  be  the  case. 

In  fact,  the  process  of  a  Taylor's  expansion  in  powers  of  one  seems  to  give 
correct  results  (on  a  transformed  space)  for  first  order  (mean  value)  results 
but  appears  to  break  down  for  second  order  results.  The  second  order  term 
2 

3  /3x3y  is  always  missing  and  seems  necessary.  This  mav  be  due  to  the  fact 
that  the  term  in  one  squared  Is  not  smaLl  compared  to  one  to  the  first  power. 

Farrell  (1976,  pp.  25,  26)  gives  a  much  improved  method  of  estimating 
Square  Law  variances  by  approximating  absorption  probabilities  and  thus  allows 
the  analysis  to  go  beyond  the  initial  stages.  He  gives  an  example  where  the 
variances  at  any  time,  t,  are  substantial. 

Clark  (1969  pp.  125,  126)  has  shown  bv  examples  that  Square  Law  variances 
start  at  zero  at  time  zero  and  tend  to  an  asymptote  at  Lfme  infinity.  In 
between  they  are  either  montonically  increasing  or.  Increasing  then 
decreasing.  Their  values  are  substantial. 

D.  The  GR  Variance 

Finally,  we  mention  that  the  only  known  CR  solution  for  more  than  one-on- 
one  (for  one-on-one,  see  reference  [1])  is  a  stochastic  duel  with  two  versus 
one  (see  Gafarlan  and  Ancker  (1984)).  A  typical  result  comparing  GR  with  SL 
is  given  in  Figure  15.  Clearly  variation  is  important,  and  GR  variance  mav 
differ  considerably  from  SL  variance. 
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Figure  15.  Two-on-One  Stochastic  Duel  with  Frlung  (2)  on  toe 
A  (Two)  Side  and  Ned  on  Lhe  b  Side 

This  section  is  concluded  with  the  observation  that  it  can  surely  be  said 
that  variance  Is  too  Important  to  be  ignored  in  any  realistic  interpretation 
of  combat  models.  This  is  a  strong  argument  for  rejecting  tine  notion  that 
some  deterministic  approximations  of  the  mean  value  of  tine  stochastic  process 
is  sufficient  in  combat  analyses. 


VI.  THE  ERRORS  IN  OTHER  MEASURES  OF1  EFFECTIVENESS 

lip  to  this  point  the  nrincipal  concern  has  been  with  Lhe  mean  value  trace 
of  the  survivors  as  a  measure  of  combat  progress.  However,  there  are  three 
other  measures  of  combat  outcome  which  are  at  least  crudely  predicted  hv  the  t, 
equations  and  which  are  now  examined,  especially  In  regard  to  their 
relationship  to  the  corresponding  SL  measures.  No  attempt  to  compare  with  the 
equivalent  GR  measures  will  be  made  here.  These  measures  are;  expected  number 
of  survivors,  expected  time-duration  of  Lhe  battle  and  the  probability  of 
winning.  These  three  predictors,  with  the  previously  discussed  mean  value 
trace  are  the  only  possible  direct  measures  that  can  he  obtained  from  the  I, 
model. 
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A.  The  Expected  amber  of  survivors 


I.  The  Linear  Law 


Weiss  (1963,  p.  598)  has  shown  that 
difference  at  Infinity  (l.e.  given  a  win 


the  conditional  SL-L  mean  value 
by  A)  is 


mA(a>) 

P(A) 


*(°°)  _ 
l 


+  b.. 


r  B  '  a0r  a _ ,Vl  _ 

la  +  6j  lat  +  8J  X  P(A) 


where  >  (8/a)b^.  Thus  L  underestimates  SL  as  this  is  always  positive. 

The  corresponding  B  side  difference  is  meaningless  since  marginal  y(°° )  =  0 
and  the  P(B)  for  L  is  also  zero  and  their  ratio  is  indeterminate.  Of  course, 
by  symmetry  the  same  can  easily  be  done  for  wins  by  B  when  <  (S/n)b^. 


However,  the  marginal  difference,  which  is  the  appropriate  one  to 
compare,  is  not  so  clear  cut.  For  example. 


again  with  a q  >  (0/a)b„.  This  will  he  very  difficult  to  explore  except  bv 
particular  examples  which  is  shown  in  Part  Two. 

2.  The  Square  and  Mixed  laws 

No  closed  form  expression  for  the  SL  mean  value  functions  at  L  =  »  exist 
although  useful  forms  for  the  marginal  distribut  Ions  of  survivors  do. 

In  general,  L  always  predicts  the  loser  with  zero  expected  survivors  and 
near  parity  this  can  be  verv  misleading.  The  discrepancy  on  the  winner's  side 
can  best  he  investigated  by  examples  as  seen  in  ParL  Two. 

It  has  been  clearly  shown  that  the  L  predictor  can  be  extremely 
misleading. 

B.  The  Expected  Tine-Duration  of  the  Battle 

Very  little  has  been  done  theoretically  on  this  measure.  However,  In 
Square  Law  parity  and  always  in  the  Linear  Law  battle  t^  +  •»  for  L.  This  is 
clearly  a  useless  result  as  all  SL  battles  have  s  finite  expected  time- 
duration  because, 

*0  b0 

f  (t)  -  l  p'(a,0,t)  +  l  p'(0,b,L) 

0  a=*0  b=0 

where  f^  =  a  proper  pdf  of  T^  (the  rv,  time-duration  of  the  combat).  The  rhs 

Is,  In  general,  a  weighted  finite  sum  of  exponentials  and  will  therefore  have 
a  finite  mean. 
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Bowen  (1965)  has  shown  that  for  any  SL  model,  the  mean  time  to  the  first 
kill  is  less  than  the  corresponding  L  time.  While  not  conclusive  concerning 
the  overall  battle,  this  is  still  an  interesting  and  suggestive  result. 

Some  particular  examples  have  been  worked  out  or  obtained  by 
simulation.  These  results  are  given  in  Part  Two  and  in  general  again  the  L 
prediction  can  be  very  misleading. 


C.  The  Probability  of  Winning 
1.  The  Linear  Law 

Brown  (1951  and  1955)  laid  the  ground  work  for  the  following  expression 
which  was  later  almost  simultaneously  puL  in  finished  form  by  Brown  (1963), 
Williams  (1963),  and  Weiss  (1963)  (most  explicitly).  For  SL, 


P(A)  =  ro/(a+B )(b0,a0)’ 


P(  B) 


where  I  is  the  Incomplete  Beta  Function  Ratio. 


'’8  /(a  +6  )  ^*0  ,b0  ^  ’ 
For  L 


P(A)  =  1, 

P(B)  =  0, 

a()  > 

(B /a )b0. 

P(A)  =  0, 

P(B )  =  1, 

a0  < 

(B/a)bn, 

P(A)  =  P(B)  =  0, 

P(D)  =  l. 

a0 

(B/a)b(1,  where  P(0)  is  n  disastrous 

draw.  The  last  expression  seems  an  appropriate  interpretation  since  both 
sides  go  to  zero  survivors  at  L  =  ”,  and  neither  ('an  be  said  to  have  won. 

From  P(A)  for  SL  we  observe  that  ()  <  T  a/(a+g)  (b^.a^)  <  1  and  that  for 

fixed  b(-j  and  a^  it  Is  monotontcally  increasing  as  n/(a+B)  increases  or  for 
fixed  a.B.b^  it  is  also  monoton  lea  1  lv  increasing  wiLh  a^.  Figure  16  Is  a 
typical  comparison  with  L. 


-37- 
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Note  that  this  general  curve  has  not  been  shown  going  through  P(A)  =*  0.5 
at  aag/8bQ  -  l  because  this  will  only  occur  for  strict  SL  parity  (i.e. 
a0  *  °o*  a  *  8),  otherwise  it  may  cross  either  above  or  below  depending  on 
!  the  exact  values  of  the  parameters.  Such  curves  are  obtained  by  fixing  ag,  bg 

and  either  a  or  8  and  varying  the  other.  It  Is  striking  that  L  can  imply  P(A) 
^  is  zero  when  it  is  >  1/2  or  P(A)  is  l  when  it  is  less  that  1/2.  The  P(B) 

|  situation  is,  of  course  symmetrical. 

'  2.  The  Square  Law 

l 

The  SL  probability  of  winning  is  also  given  in  Brown  (1951,  1955  and, 
1963)  and  Isbell  and  Marlow  (1956)  as  a  complicated  sum.  Although,  In  fact, 
it  does  have  quite  similar  properties  to  the  Linear  Law  (as  shown  hv  many 
examples  in  Part  Two)  it  is  very  difficult  to  observe  Its  properties 
theoretically. 

Brown  (1951,  1955  and  1963)  has  derived  useful  approximations  for  both 
the  SL  Linear  and  the  Square  Laws  which  are  Normal  probability  integrals, 
derived  asymptotically  but  which  are  remarkably  good  for  small  numbers.  Kisi 
,  (1966)  has  derived  the  same  expression  for  the  Square  Law  only  by  a  simple 

transformation  of  the  basic  differential-difference  equation  which  captures 
second  order  differences  (i.e.,  differences  of  differences)  and  then  replaces 
j  differences  by  derivatives. 

It  can  be  said  generally,  that  near  pari  tv  the  I.  predictor  for  winning 
can  be  extremely  misleading.  In  Part  Two  this  is  illustrated  with  particular 
examples. 


i 

i 


VII.  CONCLUSIONS 

In  the  following  we  summarize  the  main  points  in  Part  One: 

(1)  ALl  L-SL  mean  vaLue  equivalent  nairs  differ  (possiblv  considerably )  for 
all  times  except  at  crossing  points. 

(2)  At  least  for  the  Square  Law,  the  !,  trajectories  ar^  neither  a  universal 
upper  or  lower  bound  on  the  SL  mean  value  trajectories. 


(3)  Even  near  Lime  zero,  the  L  and  SL  mean  value  trajectories  mav  differ 
considerably  (they  do  not  differ  materially  for  the  Square  Law  and  the 
sequence  of  one-on-one  duels  version  of  the  Linear  Law. ) 


(4)  For  the  Linear  Law,  the  Square  Law,  the  Mfxeu  Law  and  the  Square  Law 
with  continuous  reinforcements  there  is  a  Law  of  !,arge  numbers  on  suitably 
transformed  spaces.  However  on  untransformed  spaces  one  can  only  sav 
lim  E(Alf(t)]/k  x^(t)  *  l,  Hm  E[B^(t)|/k  y  j  ( L )  =  l.  This  does  not 

k-*-oo  k-*<® 

necessarily  mean  that  as  the  Initial  force  sizes  tend  to  Infinity  the 
differences  between  L  and  SL  mean  value  trajectories  tend  to  zero.  They  mav 
even  tend  to  a  constant  or  infinity. 


(5)  Blackwell’s  Theorem  does  not  Imply  that  Individual  combatants  (and  thus 
their  superposition)  with  general  Ifts  tend  to  have  ned  Ifts  (even  after  a 
long  time).  This  Is  even  more  strongly  the  case  for  terminating  processes. 
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(6)  The  Palra-Khintchine  Theorem  does  not  imply  that  superposing  a  large 
number  of  combatants  with  general  ifts  will  yield  a  process  with  ned  ifts. 

This  can  only  be  approximately  correct  for  large  numbers  and  for  very  large 
lft  means.  Again  the  theorem  is  only  valid  for  non-terminating  processes. 

(7)  Nonhomogeneous  Poisson  processes  do  not,  in  general,  approximate  general 
renewal  processes. 

(8)  The  SL  process  variances  are  generally  quite  significant  and  can  he 
important  for  large  force  sizes,  even  near  time  zero.  In  addition,  GR  process 
variances  are  significantly  different  than  SI,  process  variances. 

(9)  The  other  L  measures,  (a)  expected  number  of  survivors,  (b)  expected 
time  duration  of  the  battle  and  (c)  probability  of  winning  are  even  less 
reliable  predictors  than  the  mean  value  trace. 

(10)  Finally,  we  emphasize  that  the  basic  -issumpt ions  of  the  SL  (and  OR  for 
that  matter)  models  simply  can  not  hoLd  for  large  numbers  of  combatants. 
Terrain  compartmentalizatlon,  weapon  ranges,  terrain  obstacles,  weather  and 
many  other  factors  (including  tactical  ones)  cause  large  scale  battles  to  he  a 
set  of  sequential  and/or  parallel  small  scale  engagements.  The  effect  of  Lhis 
point  is  illustrated  in  Figure  17,  where  one  large  battle  with  %  on  each  side 
is  compared  to  16  simultaneous  battles  of  h  on  each  side.  The  Lanchester 
solution  is  identical  for  both  cases  hut  Lhe  SL  solution  is  quite  different 


ioo(- 


Flgure  17.  A  Comparison  of  One  Large  Battle  with  Several 
Simultaneous  Smaller  Battles 


for  the  two  cases  with  the  terminal  number  of  survivors  in  the  simultaneous 
small  battles  being  80.4%  greater  than  the  one  large  battle.  Further  results 
(Ancker  and  Gafarian  (*))  are  contained  in  Section  6,  Part  Twc.  tv,,  .  5  the 
development  of  small  scale  GR  models  and  their  integration  into  large  models 
is  intrinsically  necessary.  The  effect  of  this  fractionation  on  overall 
measures  of  outcomes  is  unknown.  For  a  discussion  of  current  thinking  on  this 
matter  in  England  and  the  United  States  see  reference  (3). 


VIII.  RECOMMENDATIONS 

Both  analytical  (mathematical)  and  simulation  research  is  needed  and  is 
recommended  on  the  following  topics  on  GR  models: 

(1)  Solutions  for  moderate  size  an  and  Sq. 

(2)  Good  approximations  for  moderate  size  a^  and  bQ. 

(3)  Superposition  of  terminating  renewal  processes  (this  would  probably  be 
best  started  on  non-terminating  processes  and  proceed  towards  terminating 
processes  later). 

(4)  The  possibility  of  using  dn(t)/dt  as  the  instantaneous  rate  for  moderate 
numbers  of  superposed  iid  GR  processes. 

(5)  Numerical  techniques  to  solv<  the  complicated  analytical  models. 

(b)  In  simulations,  variance  reduction  techniques  on  non-c 1  ass  tea l 
terminating  processes. 

(7)  Determination  of  error  bounds  on  approximations. 

(8)  Integration  of  small  (or  moderate)  size  battle  models  into  large  models 
of  combat. 
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are  catalogued  under  AE  in  the  Annotated  Bibliography.  Considerable  effort 
has  been  expended  in  redrawing  all  figures  and  retyping  all  tables  in  common 
notation.  These  have  been  carefully  organized  and  then  cross  referenced 
between  the  Annotated  Bibliography  and  the  figures  and  tables  of  Part  Two  and 
between  the  sections  within  Part  Two  so  that  working  analysts  may  readily 
locate  all  the  details  they  mcy  be  interested  in. 

The  only  omitted  material  is  from  Craig  (1975)  pages  71-81,  83,  85-91, 

94,  96,  and  98-126,  and  is  a  number  of  figures  where  Square  Law  L  functions 
and  SL  mean  value  functions  are  plotted  in  the  manner  of  the  figures  in 
Section  l  following.  However,  in  the  Craig  plots  (all  for  annihilation 
situations)  all  functions  have  been  normalized  by  dividing  them  by  their 
initial  values.  This  has,  in  essence,  divided  the  SL-L  difference  by  their 
common  initial  values.  For  breakpoint  situations,  the  normalizing  factors 
were  ag  -  x(lp)  and  bg  -  y(tp).  At  first  glance  these  appear  to  be 
appropriate  transformations.  In  fact,  they  obscure  the  absolute  magnitude  of 
the  difference  (by  greatly  reducing  them).  Consequent Iv,  we  have  chosen  to 
omit  these  figures  as  otherwise  it  would  require  that  they  be  recalculated  and 
replotted  on  an  unnormalized  basis,  which  did  not  appear  to  be  worth  the 
effort. 

In  many  figures  and  tables  from  original  sources,  determinist tc  values 
are  not  given.  We  have  calculated  them  and  included  the  results  to  make 
comparisons  possible.  Any  errors  are  strictly  our  responsibility. 

A.  The  Significance  of  Parity 

The  parameters  for  stochastic  parity  (which  we  define  to  he  when  P(4)  = 
P(B))  will  not  have  the  same  values  as  they  do  for  L  parity  (except  for  strict 
parity).  Still  when  one  is  far  from  L  paritv  the  fire-fights  will  he 
decisively  lopsided  in  favor  of  the  stronger  side  in  both  models  and  not  he 
very  Interesting.  This  fact  is  especially  important  since  I.  parity  is  easv  to 
calculate  and  SL  parity  is  usually  verv  difficult  to  determine.  Strict  paritv 
is  defined  for  the  Square  and  Linear  Laws  as  a  =  3,  a,j  =  b,-,  and  Up  =  bp  and 

>  ) 

for  the  Mixed  Law  (with  A  linear  and  8  square)  is  a  =6,  u^  =  h“  ,  and  a  =  h~ 

Mon-strict  L  parity  is  as  follows:  for  the  Linear  Law,  a ( a  -  a  )  =  3(b  -  h.) 

2  2  2  2  of 

for  the  Square  Law,  a(a  -  a  )  =  3(b'  -  b  );  and  for  the  Mixed  Law  (A 

2  C  2  '  f 

Linear),  a(a0  -  af)  =  S(b()  -  bf). 

B.  Draw 

There  are  some  examples  In  what  follows  of  models  where  draws  are 
possible  outcomes.  This  matter  was  not  discussed  in  Part  One  and  deserves 
some  attention  here.  Let  us  suppose  that  each  side  has  a  breakpoint  (for 
sake  of  generality)  given  by  ap  and  bp  respectively.  Whichever  side  reaches 
it's  casualty  breakpoint  first  will  surrender  or  break  and  run.  Now,  let  ns 
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further  suppose  that  each  side  has  a  point  at  which  it's  casualties  are  such 
that  it  no  longer  has  the  capability  of  winning  but  are  not  so  large  as  to 
have  reached  it's  breakpoint.  Consequently,  it  would  disengage  if  the  other 
side  would  also  retire.  We  shall  denote  these  points  by  a  new  notation,  aR 
and  bp.  The  time-independent  joint  probability  state  space  for  stochastic 
models  of  Figure  4  is  reproduced  below  as  Figure  18  without  showing  all  the 
lattice  points  which  are  the  possible  system  states.  What  we  do  show  are  all 
the  subsets  of  states  which  have  common  possible  outcomes  in  the  breakpoint 
and  draw  situation. 


B  wins 


^ - r 

o  •  1 

B  can  win  | 

Either  con  win 

or  o  drawl 

or  a 

will  occur 1 
l 

draw  will  occur 

A  can  win  or  a 

draw  will  occur ^ 

A  wins 


A 


Figure  18.  The  Stochastic  Represcuitat  ion  of  draws 

ALl  shaded  areas  and  the  three  corners  wj  t  h  large  dots  contain  states  which 
cannot  occur. 

The  corresponding  deterministic  situation  is  shown  in  Figure  I'),  the 
phase-space  of  Figure  2(d). 


Figure  19.  The  Determin i st ir  Representation  of  Draws. 
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Any  phase-space  trace  lying  between  the  curves  shown  will  surely  end  in  a 
draw;  a  trace  below  both  will  result  in  an  A  win  and  one  above  both  will  be  a 
B  win. 

C.  Special  Models 

There  is  no  dicusssion  in  Part  One  of  two  models  which  are  in  the 
Annotated  Bibliography  and  whose  examples  are  in  this  part. 

1.  The  Weale  Special  Model 

This  model  differs  from  the  others  in  that  attrition  is  of  the 
form  a  =  a(c^  +  c0b)  and  S  =  b(c^  +  c^a),  where  Cj,  c?,  c-^,  and  c^  are 
constant  attrition  coefficients.  This  model  c ontains'as  special  cases  the 
Linear,  Square  and  Mixed  Laws. 

2.  The  Clark  LORSlfM  Model 


This  model  is  considerably  more  complicated  than  any  of  the  others. 
Briefly,  it  involves  the  following  comp  I icat ions : 

(1)  Several  different  groups  on  each  side  with  individual  group 
eharacterist ics. 

(2)  Different  kill  rates  hv  each  group  against  every  opposing  group. 

(1)  A  rate  at  which  each  group  detects  each  opposing  group. 

( 4 )  A  time  delay  from  detection  to  icon i s i t ion  for  all  in  each  side.  This  is 
called  a  shift  coefficient. 

(5)  Acquisition  autocorrelation  between  units  internal  '.a  each  group. 

Those  factors  ire  displayed  in  Tihle  r-7  for  partiml  ir  examples.  The  details 
of  this  model  are  considerably  more  involved  that  is  discussed  above,  hut  we 
shall  not  elaborate  further  here, 

I.  THE  DIFFERENCE  BETWEEN  SL  MEAN  VALUE  FUNCTIONS  AND  L  FUNCTIONS 

All  figures  and  tables  in  this  section  lispl  iv  the  ditference  except 
Tables  1-2,  i-l(a),  I-b  and  Figure  1-18,  whicli  . I . >  not  have  the  deterministic 

information.  They  ire  included  for  comp  I eLoness. 

Additional  information  on  the  differences  is  given  in  Figures  VI-1 
through  VI-12. 

Parity  may  be  observed  in  Figures  I-I(b),  t-1  (solid  curves  on  1 v ) , 

I-5( c  ) , ( d ) ,  I  —  21(b),  [-24(c),  I25(a),(c),(e),  l-2h(c),(d),  vr-1  and  VI-4  while 
strict  parity  occurs  in  Figures  I-l(a),(c),  1-4,  l-1),  and  1-12.  There  are 
scattered  parity  points  in  various  tables  hut  being  at  isolated  Lime  points 
they  are  of  Little  interest.  As  observed  earlier,  at  or  near  parity  L  is 
particularly  poor  as  an  approximation  to  SL. 
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An  important  point  about  breakpoints  is  illustrated  in  Figures  1-3,  1-4, 
and  1-5.  That  is,  for  L  the  breakpoint  curves  are  identical  with  the  non¬ 
breakpoint  curves  until  the  breakpoints  are  reached,  at  which  point  the  L 
curves  are  horizontal  straight  lines.  However,  the  SL  breakpoint  curves  are 
higher  than  their  non-breakpoint  counterparts  and  depart  from  the  L  curves 
much  earlier  in  time.  Contrary  to  some  speculation,  breakpoints  do  not 
improve  the  L  approximation  to  SL  mean  value  functions;  the  effect  of 
absorption  probabilities  simply  occurs  earlier. 

In  Figures  1-19  through  1-26  we  see  a  two  versus  one  model  illustrate  the 
differences  in  L,  SL,  and  GR  mean  value  functions.  One  striking  point  is  that 
for  the  same  value  of  a  =  (  p^  /  U  various  values  of  greatly  change  the 
GR  curve.  In  other  words,  combining  the  parameters  and  can  be  grossly 
inadequate. 

The  other  points  mentioned  on  pages  15  and  lo  are  well  illustrated  by 
these  results,  and  they  emphasize  that: 

(1)  L  functions  are  generally  Inadequate  as  an  approximation  to  SL  mean 
value  functions,  especially  in  the  most  interesting  cases  where  one  side  does 
not  have  a  lopsided  preponderance  of  force. 

(2)  SL  mean  value  functions  are  generally  inadequate  as  an  approximation  to 
their  GR  equivalents  (at  least  for  the  small  fire-fights  so  far  considered). 


12 


JAMES  (1981)  SQUARE  LAW  p.S2 


TIUE 


JAMES  (1981)  SQUARE  LAW  p.52 

aQ  =12,  b0  =  8 


Figure  1-3 
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a  SURVIVORS 


Figure  1-8. 
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GYE  AND  LEWIS  (  1976)  p. 118 
SQUARE  LAW 


a0  =  23,  b0  =  32 

af  *  bf  =  0 
a  =  8 


TIME 

'"A(t).  "iBft) 

L  POINT 
WITH  SAME  X 

TIME 

TAl't).  iRB(t) 

L  POINT 
WITH  SAME  X 

3.1 

22,  32 

22,  31.3 

46.9 

10,  25 

10,  24.4 

9.7 

20,  30 

20,  29.9 

55.2 

8,  24 

3,  23.6 

16.5 

18,  29 

13,  23. 6 

63.9 

6,  24 

6,  23.0 

23.7 

16,  28 

16.  27.4 

72.6 

4.  23 

4,  22.6 

31.1 

14,  27 

14,  26.3 

31.6 

2,  23 

2,  22.3 

38.9 

12,  26 

12,  25.3 

CLARK  ' 1 969)  p.!22 
SQUARE  LAW 


*f  »  bf  »  0 


*0 

ho 

a 

8 

t 

■A(t> 

^l(t) 

*(t) 

y(t) 

AA(t> 

AB(t) 

6 

— 

6 

.001 

.004 

275. 

1.059 

5.134 

0.000 

5.195 

1.059 

-.061 

6 

6 

.0015 

.004 

291. 

1.133 

4.640 

0.000 

4.740 

1.133 

-.100 

6 

6 

.002 

.004 

312. 

1.227 

4.102 

0.000 

4.*M0 

1.227 
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6 

6 

.004 

.004 

2500. 

1.929 

1  .929 

0.000 

0.000 

1.929 

1.929 

8 

a 

.001 

.004 

275. 

1  .226 

6.363 

0.000 

6.922 

1 .226 

-.059 

8 

8 

.0015 

.004 

291. 

1.314 

6.213 

0.000 

6.320 

1.314 

-.107 

8 

8 

.002 

.004 

312. 

1.426 

5.492 

0.000 

5  623 

1.426 

-.136 

3 

8 

.004 

.004 

2500. 

2.407 

2.407 

0.000 

0.000 

2.407 

2.407 

12 

6 

.001 

.004 

2500. 

3.297 

1.969 

0.081 

0.040 

3.216 

1.929 

12 

6 

.0015 

.004 

493. 

6.474 

1.339 

6.920 

0.000 

-.446 

1.339 

12 

6 

.002 

.004 

312. 

8.160 

1.176 

8.480 

0.000 

-.320 

1.176 

12 

6 

.004 

.004 

137. 

10.275 

1.060 

10.390 

0.000 

-.115 

1.060 

TabLe  1-1 


CLARK  ( 1969)  p. 153 

SQUARE  LAM  'WITH  STOCHASTIC  ACQUISITION) 


af  =  bf  =>  0 


m 

Ml 

'Bui 

n 

0 

L 
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x(t) 

y(t) 
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1 

5 
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.312 

.006 

5 

6 

.00A 

.004 
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\ 
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.004 
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.013 
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6.383 
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! 
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0.000 

4.568 
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) 
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.004 
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.004 
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’  1 

6 
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1  .504 

.408 

.608 
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1 2 
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.  ?'*•* 

3.364 

.790 

.785 

1 2 
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8.268 
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.061 

Astj  isieion  Probabilities  -  (p  =  <7  =  .85). 


Table  l- A. 
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KARR  (1975b)  LINEAR  LAW  p.10 


Figure  1-11. 
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"A(t) 
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0.004 
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SPRINGALL  (1968)  MODEL  p.164 

«qs40,  bQ  =30,  m j  =2  9,  r»2  =20,  O  =0. 1  0, 
flsO.1,  y  =0.2,  JsO.  1 


—  L 
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SPRINGALL  (1968)  MODEL  p.165 


b q  =4  0,  m2  =50,  na=40,  Q si,  8  =  1,  y  =10,  i=10 
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Figure  1-13. 
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CLARK  (1982)  pp. 157, 158 
LORSUM  MOOEL 


This  table  defines  the  parameters  for  following 
Figures  1-14  through  1-17  and  Table  1-8. 


CASE 

INPUT  OATA 

. 

1 

Initial  A  Group  1  Weapons 

2 

Initial  A  Group  2  Weapons 

10 

Initial  8  Group  1  Weapons 

1 

Initial  B  Group  2  Weapons 

5 

i 

Rate  A  Group  1  Weapons  Kill  Acquired 

8 

Group  1 

Weapons 

.001  : 

Rate  A  Group  1  Weapons  Kill  Acquired 

B 

Group  2 

Weapons 

.002 

Rate  A  Group  2  Weapons  Kill  Acquired 

8 

Group  1 

Weapons 

.0002  ; 

; 

Rate  A  Group  2  Weapons  Kill  Acquired 

8 

Group  2  Weapons 

.0004 

Rate  8  Group  1  Weapons  Kill  Acquired 

A 

Group  1 

Weapons 

.004  ! 

Rate  8  Group  1  Weapons  Kill  Acquired 

A 

Group  2 

Weapons 

.008  ■ 

Rate  8  Group  2  Weapons  Kill  Acquired 

A 

Group  1 

Weapons 

.0008 

' 

Rate  8  Group  2  Weapons  Kill  Acquired 

A 

Group  2 

Weapons 

.0016  ■ 

A  Shi  ft  Coefficient 

10.0 

8  Shift  Coefficient 

12.0 

Rate  an  A  Detects  a  Firing  8 

.03 

Rate  a  8  Detects  a  Firing  A 

.05 

Rate  an  A  Detects  a  Silent  8 

.015  ; 

Rate  a  8  Detects  a  Silent  A 

.025 

Rate  an  A  Loses  a  Detected  8 

.01 

Rate  a  8  Loses  a  Detected  A 

.01 

A  Observer  Autocorrelation 

0.0 

8  Observer  Autocorrelation 

0.0 

2 

Initial  A's  are  Reduced  by  50?  in  Each 

Group 

3 

Each  A  Kill  Rate  is  Reduced  by  25? 

4 

Rate  A  Oetects  a  Firing  B 

.005 

Rate  B  Oetects  a  Firing  A 

.005 

Rate  A  Oetects  a  Silent  B 

.004  , 

; 

Rate  8  Oetects  a  Silent  A 

.004  j 

5 

Observer  Autocorrelations 

1 

.25  ; 

t 

6 

Detection  Rates  Equal  to  Case  4,  Observer  Autocorrelations 

.5 

IB  -  6B: 

Same  as  l  -  6  with  all  •jeoup  I’s  eliminated . 

Table  1-7 
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CLARK  (1982)  p.10O 


Time 

CLARK  (1982)  p.181 


Figure  I- 14 
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CLARK  (1982)  p.l«4 


CLARK  (1982)  p.  165 


Figure  I-L6 


CLARK  (1982)  p.168 
LORSUM  MODEL 


CASE  2 


CLARK  (1982)  p.1«7 
LORSUM  MODEL 
CASE  2 


00  200.0  400.0  000.0  000.0  1000.0 

TIMO 


(a) 


(b) 


Figure  1-17 
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CLARK  (1982)  pp. 168-169 
LORSUM  MODEL 


t  =  500 


A  GROUP 

1  SURVIVORS 

B  GROUP 

1  SURVIVORS 

DETERMINISTIC 

DETERMINISTIC 

CASE 

STOCHASTIC 

ERROR 

STOCHASTIC 

ERROR 

.5046 

-.0828 

.4712 

.2290 

.1126 

-.0879 

.7659 

.1298 

.4311 

-.2494 

.5678 

.2184 

1.0179 

.0351 

.7018 

.0977 

5 

.5274 

-.2515 

.4689 

.2290 

6 

1.1318 

-.0011 

.6932 

.1011 

A  GROUP 

2  SURVIVORS 

B  GROUP 

2  SURVIVORS 

DETERMINISTIC 

DETERMINISTIC 

CASE 

STOCHASTIC 

ERROR 

STOCHASTIC 

ERROR 

1 

7.217 

.602 

3.535 

-2.757 

2 

1.202 

-  .654 

4.525 

-  .850 

3 

6.948 

.626 

3.979 

-2.083 

4 

6.292 

.425 

3.463 

-2.366 

5 

7.183 

.814 

3.571 

-2.573 

6 

6.824 

.247 

3.787 

-2.092 

1H 

6.807 

2.152 

3.492 

-3.490 

2H 

1.757 

1.084 

4.379 

-3.309 

3H 

6.535 

2.080 

3.827 

-3.792 

4H 

6.960 

1.507 

3.837 

-3.402 

5H 

6.694 

2.216 

3.489 

-3.487 

6H 

7.499 

1.014 

4.103 

-3.283 

m  m 


Figure  1-22 


Figure  L-24 
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00  ?  00 


Figure  1-25. 
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GAFARIAN  and  ANCKER  (1984)  SQUARE  LAW  pJ-10  GAFARIAN  and  ANCKER  (1984)  SQUARE  LAW  pJrll 


Figure  1-26. 
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II.  THE  EXPECTED  NUMBER  OF  TERMINAL  SURVIVORS 


This  section  is  closely  related  to  Section  I  and  in  fact  is  exactly  the 
situation  examined  there  at  t-MB  only.  It  Is  singled  out  for  particular 
attention  because  significant  discrepancies  between  terminal  L  and  SL  (or  GR) 
mean  values  should  be  especially  important  to  decision-makers  and  other  users. 

Most  of  the  figures  and  tables  in  Section  I  and  Figures  VI-3  through  VI- 
12  have  been  carried  far  enough  in  time  so  that  terminal  differences  are  very 
closely  approximated.  This  is  true  except  for  the  following: 

(L)  Where  there  is  no  deteministic  information  given;  i.e.  Figure  1-8  and 
Tables  1-2,  I-3(a)  and  1-6. 

(2)  Where  time  is  not  carried  far  enough;  i.e.  Figures  1-9  and  1-12  and 
Table  I-3(b). 

In  addition  to  the  parity  cases  noted  in  Section  I  we  also  have  parity 
here  in  Figures  II-2  and  II-3.  Again  we  note  that  x(°°)  =  y(»)  =  0. 

The  Inadequacies  of  L  and  SL  models  as  previously  noted  on  pages  15  and 
16  and  in  Part  Two,  Section  I  are  even  more  apparent  here. 
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Or'fi'  |«.7f’  WALL  I 


WEISS  (1963)  LINEAR  LAW  p.601 


Figure  II-4 


KISI  (1965)  p.52 
.  LINEAR  LAU 

*0  *  3*  bo  =  2 

af  *  bf  =  0 
a  =  0 


LAU 

u 

P(A) 

E[A(-)|A] 

SL 

Linear 

1  1  3 

4  4  16 

0.688 

2.09 

L  Linear 

1.0 

.  1.0 

Table  II-3 
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SPRINT.AU  (I96H)  P.  158 
SPRINGALL  NODEL 


•»  •  5,  a  »  0.9,  B  *  3.0.  a  *  Z.o,  y  •  2.0 

*r  *  br  *  0 

t  •  - 


H 

«•) 

EXPECTED  OMOI  OF  B  SUttIWSS,  b 

expected  meat  of  a  sum*c«s,  < 

OCTOVHRISnC 

STOCHASTIC 

OCTERMINISTIC 

STOCHASTIC 

5 

0.88507 

2.500 

2.842 

0 

0.203 

10 

0.97080 

5.S87 

S.854 

0 

0.059 

zo 

0.99763 

11.762 

12.003 

0 

0.005 

30 

0.99979 

17.937 

18.117 

0 

0.001 

40 

0.99998 

24.112 

24.263* 

0 

0.000 

50 

1.00000 

30.287 

30.529 

0 

0.000 

60 

1.00000 

tt.462 

36.706 

0 

0.000 

SPRINGALl  II96H)  p.159 
SPR INGALL  HOPEI 


«q  -  50,  IK,  -  40,  •  50,  a  •  fl  •  V,  y  %  5  -  10 


IB 

EXPECTED  AMEX  OF  SUttlVCXS  «„(-) 

DCTEfttlMSTlC 

STOCHASTIC 

0 

0.14455 

0 

0.704 

I 

0.14180 

I 

1.680 

2 

0.U9OO 

2.656 

4 

0.13321 

4 

4.608 

6 

0.12715 

6 

6.560 

3 

0.12082 

3 

8.512 

10 

0.U42O 

10 

10.465 

20 

0.07620 

20 

20.243 

JO 

0.03071 

30 

30.066 

SPRINGALL  '  1  968  !  p.177 
SPRINGALL  HQOCL 


*o  •  60.  h,  •  60.  a 


0.9, 

4f  ‘ 


6  •  1.0,  y  •  2.0, 

bf  •  5 


2.0 


ED3 

•»<*> 

V> 

m 

ca 

10 

7.31188 

11 .68936 

0 

9.64 

15 

7.2242? 

11.84205 

D 

9.85 

20 

7.18098 

11.91717 

m 

9.95 

30 

7.14259 

11.98471 

10.04 

40 

7.17821 

i?.mni6 

B 

10.07 

50 

7.1228? 

17.01975 

B 

10.08 

60 

7.12149 

l?.0??I2 

H 

10.08 

(b) 


(c) 


Table  1 1 -4 
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III.  THE  EXPECTED  TDC-PORATIOH  OF  THE  COMBAT 


The  following  do  not  give  L  termination  times  for  comparison  with  E(TQ] 
but  are  included  for  completeness;  Tables  III-2  and  III-3. 

Certain  figures  in  Section  I  contain  information  on  time-duration  and 
should  be  consulted  in  connection  with  the  material  given  here.  They  are 
Figures  I-L(b),  (c)  (these  are  parity  cases  where  t^  ®),  1-2,  1-6,  1-7,  and 

1-8  which  are  all  Square  Law. 

The  following  figures  clearly  indicate  the  effect  of  breakpoints  when 
compared  to  annihilation  models:  Figures  IEI-L(b),  tll-2(b),  lll-3(a),  [II-4 
and  Table  Ill-l(b).  The  breakpoints  lower  the  time-duration  of  the  fire-fight 
substantially  as  one  would  expect. 

Nowhere  is -the  variance  of  Tp  explicitly  computed  except  in  Tables  TIC  — 
1(b),  (c).  The  5th  and  95th  percentiles  are  shown  in  Table  III— 1 (b). 

However,  these  tables  and  an  inspection  of  all  the  probability  distribution 
functions  (fairly  flat)  and  density  functions  (rather  spread  out)  demonstrate 
that  v[Tq]  is  indeed  substantial.  This  fact  and  the  large  discrepancies 
between  b.[Tp]  and  tp  which  can  he  extremely  large  near  parity  (where 
tf  +  «)  lead  to  the  conclusion  that  tp  is  an  extremely  poor  approximation  to 
tne  measure,  expected  time-duration. 

Finally,  Figures  1 1 1  —  1  and  TIC  —2  not  only  show  plots  of  the  distribution 
function  of  TR  but  also  the  marginal  functions  P[TR  <  t,  A  wins]  and  P[Tp  <  t, 
B  wins].  We  point  out  that  these  latter  two  are  not  conditionals  (some  of 
which  are  shown  in  Tables  I I X— 2  and  III— 3(a)),  but  are  improper  distribution 
functions  with  the  property  that  as  t  *  ®  these  function  Lend  to 

*0  b0 

1  p(a,0,*)  =  P(A)  and  £  p  ( 0 ,  b ,°° )  =  P(8).  This  is  related  to  the  fact  that 

a=ap  b=bp 

the  expression  at  the  bottom  of  page  36  is  the  time  derivative  of  the 
identity; 

Ft  (t)  =  P[Tq  <  t]  =  P(T0  <  t,  A  wins]  +  PlT^  <  t,  B  wins] 

ao  bn 

=  l  p(a,0,t)  +  l  p( 0 , b , t ) 

a“af  h=bf 

which,  of  course,  only  involves  the  absorption  probabilities. 
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JAMES  (1981)  SQUARE  LAW  p.53 


2  a  4  5  8  7  S  9  10  11  12 

TIME 


JAMES  (1981)  SQUARE  LAW  p.54 


Figure  Ill-l 


JAMES  (1981)  SQUARE  LAW  p.24 


TIME 


JAMES  (1981)  SQUARE  LAW  p.25 


*o  =10 
bo  =10 

-f=bf=7 
a  =  $ — 0.05 


90% 


WEALE  (1976)  pp. 43-45 
SQUARE  LAW  (WITH  DRAWS) 


ao  =  b0  =  10 

af  ■  2,  bf  •  3 

a  =  0.16,  8  =  0.08 

P(A)  =  0.77630301,  P(B)  *  0.03090876,  P(D)  *  0.19278822 

DRAWS  OCCUR  IF: 

a  *  2,  b  <  6  or 
b  »  3,  a  <  5 


E[T0]  = 

5.35 

ECTD|a] 

=  4.58 

E[T0|b]  = 

6.21 

e[Td|d]  ■ 

8.31 

TJ 

ft0 

fT0 

Ftd|a 

m 

FTd|B 

flD|B 

FTd|D 

f  tD  |  D 

0.00000 

0.00000 

0.00000 

0.000 

0.00000 

2.00 

0.03466 

0.073 

0.04445 

0.093 

0.00419 

0.012 

0.00010 

4.00 

0.33787 

0.193 

0.42414 

0.234 

0.14074 

0.140 

0.02210 

0.035 

0.67320 

0.130 

0.79953 

0.127 

0.50572 

0.192 

0.19136 

0.134 

8.00 

0.86077 

0.063 

0.95205 

0.038 

0.81202 

0.107 

0.50103 

0.157 

10.00 

0.94620 

0.027 

0.99118 

0.008 

0.94756 

0.037 

0.76484 

0.102 

12.00 

0.98115 

0.010 

0.99866 

0.001 

0.98824 

0.009 

0.90952 

0.047 

14.00 

0.99400 

0.004 

0.99982 

0.000 

0.99774 

0.002 

0.96995 

0.017 

16.00 

0.99825 

0.001 

0.99998 

0.000 

0.99961 

0.000 

0.99104 

0.006 

PP1 

0.99952 

0.000 

1.00000 

0.000 

0.99994 

0.000 

0.99754 

0.002 

Wmm 

0.99988 

0.000 

1.00000 

0.000 

0.99999 

0.000 

0.99936 

0.000 

mm, 

0.99997 

0.000 

1.00000 

0.000 

1.00000 

0.000 

0.99984 

0.000 

24.00 

0.99999 

0.000 

1.00000 

0.000 

1.00000 

0.000 

0.99996 

0.000 

26.00 

1.00000 

0.000 

1.00000 

0.000 

1.00000 

0.000 

0.99999 

0.000 

28.00 

1.00000 

0.000 

1.00000 

0.000 

1.00000 

0.000 

1.00000 

0.000 

30.00 

1.00000 

0.000 

1.00000 

0.000 

1.00000 

0.000 

1.00000 

0.000 

Table  III-2. 
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SPnmGALL  (1968)  MODEL  p.161 


,„.o.  a  =o..,y -i.o.  S=2.o, 


nasbo 


SPRINGALL  (1968)  MODEL  P-162 


•  o  or  b0 

Figure 


WEALE  (1976)  oo.15.I6.t7 

SPECIAL  MODEL 


*,  •  b8  •  10 
*r  *  "f  *  3 

P(A)  -  .99700850,  P(8)  •  .00001828,  P(0)  •  .00717318 

DRAMS  OCCUR  IF: 

*  ■  2.  b  *  6  or 
b  «  3.  »  <  5 


AIT1ITI0R  FURCTIORS: 

A  Sid*  •  *10.05  ►  0.05b) 

8  Sfd*  •  b(0.02S  *  0.005*) 


1  ECrj] 

2.02 

£tT0|*J 

•  2.00 

ECToIrJ  •  3.43 

ElToloJ 

4.44 

T 

fT0 

Jr° 

FT0|A 

fTo|A 

FTo|8 

fTO|0 

Ft0|D 

fTD|D 

a.oo 

Q.QOOOQ 

0.000 

^  1 

0.00000 

0.000 

0.00000 

0.000  ‘ 

0.00000 ' 

0.000 

1.00 

0.08261 

0.314 

0.00320 

0.316 

0.00235 

0.015 

0.00009 

0.001 

2.00 

0.56262 

0.471 

0.56656 

0.474 

0.09754 

0.210 

0.01865 

0.059 

3.00 

0.87641 

0.172 

0.80158 

0.171 

0.40279 

0.353 

0. 16248 

0.231 

4.00 

0.97035 

0.043 

0.97420 

0.041 

0.71639 

0.250 

0.43818 

0.290 

5.00 

0.99267 

O.QIO 

0.99483 

0.008 

0.89422 

0.114 

0.69471 

0.211 

6.00 

0.99797 

0.002 

0.99900 

0.002 

0.96638 

0.041 

0.85568 

0.115 

7.00 

0.99936 

0.001 

0.99981 

0.000 

0.99038 

a.  012 

0.93695 

0.054 

8.00 

0.99977 

0.000 

0.99996 

0.000 

0.99743 

0.003 

0.973S0 

0.023 

9.00 

0.99991 

a.ooo 

0.99999 

0.000 

0.99934 

0.001 

0.98904 

0.010 

10.00 

0.99997 

0.000 

1.00000 

0.000 

0.9998* 

0.000 

0.99549 

a. oo4 

11.00 

0.99999 

0.000 

1.00000 

0.000 

0.99996 

0.000 

0.99814 

0.002 

12.00 

0.99999 

0.000 

1.00000 

0.000 

0.99999 

0.000 

0.99924 

o.ooi 

13.00 

1.00000 

0.000 

1.00000 

0.000 

1.00000 

0.000 

0.99969 

0.000 

14.00 

l. 00000 

0.000 

i.ooooo 

0.000 

1.00000 

0.000 

0.99987 

0.000 

15.00 

4.00000 

0.000 

1.00000 

0.000 

1.00000 

0.000 

0.99995 

0.000 

16.00 

1.00000 

0.000 

1.00000 

0.000 

1.00000 

0.000 

0.99990 

0.000 

SPPrw.AU  (1968)  ».  177 
5P9I9GALL  9Q0EL 


*o  •  60.  bg  *  60,  o  ■  0.9,  g  •  1.0,  T  *  2.0,  $  -  2.0 
*f  •  bf  •  5 


»2  4*6 

Efr0J 

10 

0.45087 

15 

0.23260 

20 

0.15671 

30 

0.10796 

40 

0.09529 

50 

0.09165 

60 

0.09091 

Table  III-3. 
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IV.  THE  PROBABILITY  Of  WINNING 


We  note  again  that,  for  L  parity,  the  L  prediction  Is  that  both  sides  are 
annihilated  or  both  sides  reach  their  breakpoints.  In  either  case,  neither 
side  wins  and  therefore  L  predicts  P(A)  =  P(8)  =  0.  At  parity,  the  L 
probabilities  have  a  discontinuity  and  make  a  jump  with  a  predicted 
probability  of  zero  at  the  parity  point.  At  the  corresponding  SL  point,  the 
win  probabilities  are  P(A)  =  P(B)  *  i/2  for  strict  L  parity  and  will  be  some 
value  different  from  1/2  for  non-strlct  L  parity.  In  the  latter  case,  it  may 
even  be  on  the  wrong  side  of  1/2,  e.g.,  SL  P(A)  may  be  greater  than  1/2  when  L 
P(A)  =*  0.  The  jump  points  are  shown  as  dashed  vertical  lines  on  almost  all 
curves  or  are  otherwise  indicated  on  some  curves  and  all  tables  except  Figure 
IV-15  and  Table  IV-10  where  no  L  information  is  given. 

Some  other  sections  in  Part  Two  contain  figures  and  tables  on  the 
probability  of  winning  and  should  be  consulted.  These  are: 

(1)  The  Square  Law;  Figures  I-l(b),  (c),  1-2,  1-6,  1-7,  1-8,  tll-l,  1II-2 
(in  these  latter  two  figures  F  (« )  are  P(A),  P(B)  on  the  A,  B  victory  curves 

D 

respectively),  V-12,  V-13,  V-14  and  Tables  II-l,  II-2,  and  V-9. 

(2)  The  Square  Law  with  draws;  Table  111-2 

(3)  The  Linear  Law;  Tables  l-5(a)  and  1 1  —  3 . 

(4)  The  Springall  Model;  Tables  ll-4(c),  V- 6  and  V— 11. 

(5)  The  Weale  Special  Model;  Table  llt-3(a)  and  V-10. 

Some  of  these  contain  L  comparisions  and  some  do  not.  Usually,  the  L 
predictions  are  obvious  or  easily  determined  from  the  discussion  in  Part  One. 

Tlie  probability  of  winning  figures  all  are  shown  as  continuous  functions 
of  the  parameters  ag,  bg  (also  ap  and  bp  in  some  cases),  a,  S  (which  may  be 
exhibited  in  terms  of  PA,  Pfi,  u^,  and  u^).  If  af)  nnd  bg  are  fiKer|>  anfi 
if  a  (or  3)  or  some  function  of  a  and  B  is  varied  the  curves  are  genuinely 
continuous.  However,  this  is  the  case  only  for  Figures  IV— l l ,  IV-10(b) 
(showing  slope  discontinuities  which  it  should  not)  and  IV— 17.  All  others 
should  have  values  at  certain  discontinuous  points  only.  In  other  words,  if 
ag  or  bg  is  the  variable,  the  P(A)  and  P(B)  curves  should  he  discontinuous. 
This  can  be  somewhat  misleading  and  care  should  be  taken  in  reading  the 
curves.  Figure  IV-2  illustrates  this  point,  a  and  B  are  fixed,  ap  =  bp  =  0 
and  the  curve  parameters  are  ag  +  bg.  If,  for  example,  ag  +  bg  =  o,  ag  can 
onLy  assume  values  0,1, 2, 3, 4,  and  5  while  the  corresponding  values  of  bg  are 
5, 4, 3, 2, l,  and  0.  P(8)  only  has  non-zero  values  at  abscissa  points  0,  1/16, 

4/9,  9/4,  4  and  130  .  Interestingly,  there  is  only  one  value  to  the  right  of 
the  jump  point.  Also  one  has  to  be  curious  about  the  curve  whose  parameter  is 
a0  +  bQ  *  0  as  this  implies  neither  side  has  any  combatants!  If  the  authors 
used  some  limiting  process  to  obtain  this  curve  or  assumed  ag,  bg  continuous 
for  the  others,  they  do  not  say  so. 

In  any  event,  the  presented  material  contains  a  wealth  of  evidence  that  L 
is  a  very  poor  predictor  of  the  SL  probability  of  winning. 
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.EE 


LEE  and  WANNASILPA  (1972)  p.29 


LEE  and  WANNASILPA  (1972)  p.30 


P(B) 


SQUARE  LAW: 


Figure  TV- 5. 
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CRAIG  (1975)  SQUARE  LAW  p.47 


£  <■», 


KARR  (1975a)  SQUARE  LAW  p.18 


DATHE  (1967)  SQUARE  LAW  p.  1  8  (FIG.  2) 


Figure  IV—  1 2 


-119- 


BROKEN  VERTICAL  LINES  ARE  JUMP  POINTS  FOR  L  P(A) 


DATHE  ( 1967)  SQUARE  LAW 

p.20  Fig  4 


N  (.  0 A  BHEO  VERTICAL  LINES  ARE 
MODEL  JUMP  POINT!  PROM  P|AI  =  0  TO  I 


DATHE  (1967)  SQUARE  LAW  p.21,  Fig  5 


LEE  (1979)  p. 10 

SQUARE  LAW 


1. 

0 

2. 

3. 

6.0 

7.0 

1.0 

2.0 

3.0 

4.0 

5.0 

2.0 

1.0 

7.0 

6.0 

5.0 

4.0 

3.0 

6.0 

7.0 

1.0 

1.0 

7.0 

2.0 

6.0 

3.0 

4.0 

5.0 

4.0 

5.0 

3.0 

6.0 

2.0 

7.0 

1.0 

0000248 


0081507 


1126240 


5000000 


3873759 


9938492 


9999751 


0007054 

.0003894 


.0187682 


.8083533 

.9812317 

.9996105 

.9992945 

.0036160 

.0000602 

.0050362 

.9146103 

.9949637 

.9999397 

.9963839 


KARR  (  1976)  p..’2 

SQUARE  LAW 


if  *  bf  ■  0 


*  -  50;  b  •  50 

«  •  200;  b  -  200 

.77 

.91 

.74 

.79 

.54 

.68 

.42 

.34 

.23 

.25 

.20 

.08 

T.ihle  IV- 2 


GTE  AMO  LEWIS  (  1976)  p.  114 
SQUARE  LAW 


LEE  AMO  HANNAS ILRA  (197?)  0.25 
SQUARE  LAW 


Af  *  bf  *  0 

a  1 
8  *  J5 


if  *  bf  •  0 


INITIAL 

FORCES 

0(8) 

*o 

Op 

»  •  .OS.  •  -  .05 

S  »  .1,  •  •  .3 

9  *  .05,  a  ■  .2 

1 

3 

0.9583 

- 

- 

2 

6 

0.9938 

0.37179 

0.64338 

3 

9 

0.9990 

0.92326 

0.84266 

4 

12 

0.9998 

0.95268 

0.38234 

s 

IS 

0.9999 

0.97028 

0.91067 

a 

18 

1.0 

0.98110 

0.93148 

7 

21 

1.0 

0.98787 

0.94703 

8 

24 

1.0 

0.99219 

0.95881 

9 

27 

1.0 

0.99554 

0.96804 

10 

30 

1.0 

0.98188 

20 

60 

1.0 

* 

• 

L  for  0(8) 

1.0 

l  .0 

1.0 

(a) 


P(B) 

*0 

HHB 

2 

3 

4 

5 

i 

U.962 

0.999 

7 

0.390 

0.995 

1.000 

1.000 

3 

0.795 

0.984 

0.  '99 

1.000 

4 

0.6U5 

0.962 

0.997 

.  .000 

s 

U.S71 

0.926 

0.993 

1.000 

6 

0.461 

o.a/6 

0.984 

0.999 

7 

0.  360 

0.313 

0.970 

0.997 

3 

0.233 

0.947 

0.993 

9 

0.2'M) 

0.656 

0.916 

0.937 

10 

0.143 

0.571 

0.075 

0.977 

11 

0.099 

0.486 

0.326 

0.962 

12 

0.i)67 

0.405 

0.768 

0.941 

13 

0.044 

0.330 

0.703 

0.414 

14 

0.020 

•0.264 

0.634 

0.3/9 

15 

0.018 

0.207 

0.563 

0.333 

La 

9.011 

0. 160 

0.492 

0.  790 

17 

0.006 

0.121 

0.423 

0.737 

Id 

0.004 

0.090 

0.359 

0.679 

W 

0.002 

0.066 

0.619 

20 

o.ooi 

0.047 

0.246 

0.557 

21 

0.001 

0.033 

0.200 

0.495 

j  2 

0.000 

0.023 

0.160 

0.434 

2  b 

- 

- 

- 

0.496 

mrrz . 

Above  the  he<*vi4  lines  ( f 

for  L ,  P(B)  * 

l.Ol  *nd  below  c he  lines,  P(B)  Cor  L  > 

0. 

(b) 


Table  IV-3 
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KARR  '  1 97‘3<i '  p.?l 
SQUARE  LAW 


HARK  '197*^'  p.  1  7 


For  ,  P:d 1  =  0  ab^ve 


EE 


Figure  IV-18 
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LEE  ANO  WANNASILPA  (1972)  p. 25 
LINEAR  LAW 


af  =  bf  =  0 


INITIAL 

FORCES 

P(B) 

ao 

El 

8  *  .05,  a  =  .05 

bred 

8  =  .05,  «  =  .2 

1 

3 

0.8750 

0.57812 

0.51200 

2 

6 

0.9375 

0.55505 

0.57680 

3 

9 

0.9672 

0.54480 

0.61740 

4 

12 

0.9824 

0.53871 

0.64816 

5 

15 

0.9904 

0.53458 

0.67329 

6 

18 

0.9947 

0.53158 

0.69469 

7 

21 

0.9970 

0.52917 

0.71339 

8 

24 

0.9983 

0.52727 

0.73003 

9 

27 

0.9991 

0.52570 

0.74501 

10 

30 

0 

0.9995 

0.52437 

• 

0.75864 

• 

20 

• 

60 

1.0 

• 

0.51720 

• 

0.85079 

L  for  P(B) 

1.0 

0 

0  ; 

J 

Table  IV-8 
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SPR INGALL  (1960)  pp, 168, 169 
3PR INGALL  MODEL 

5,  a  ■  0.9,  0  ■  1.0,  S  •  r  -  2.0.  if  *  t>f  ■  0 


SPRINGALL  (1968)  MODEL  p.156 


AHO  bo 

2(B) 

5 

0.53560 

10 

0.55793 

:o 

0.58710 

30 

0.60847 

40 

0.62600 

50 

0.64112 

60 

0.65453 

"2  ’  5.  „  ■  0.9,  0  ■  3.0.  s  •  T  ■  2.0,  If  *  bf  -  0 


b  *20,  m  iS,  n  «5,  •  i*0, 
0  2  2  f 

bf so,  aso.«,  8; i.o.  y  .2.0 
i  =  2.0 


20  '0  *.0 

«  - — 

0 


Figure  IV-22 


IlgUOb, 

P(B) 

1  5 

0.88567 

10 

0.97080 

:o 

0.99763 

30 

0.99979 

40 

0.99998 

c’v.-v. 

;  ',kJ  r 

■^0  !N0AL.  '*0D£, 

60.  i  0.9. 

8  '  1.0. 

V  ■ 

bf  ■  5 

!~<n,  m<1  n„ 

P(A)  - 
*-  i 

!  1  3 

1.12519 

•  L 

).  31638 

1 

3.31211 

1  >•' 

0.30832 

j  30 

0.30690 

SO 

0.30637 

I  60 

1 

0.30624 

Table  IV- 10 
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V.  THE  VARIANCE  OP  THE  PROCESS 


In  this  section  standard  deviations  are  shown  directly  in  Figures  V-l 
through  V-9  and  Tables  V-l  through  V-b.  The  remainder  of  the  section  shows 
state  probability  density  functions  and  probability  mass  functions  in  various 
forms.  The  pdf's  and  pmf's  are  given  here  even  though  they  do  not  explicitly 
display  variance  (or  standard  deviation)  because  they  visually  illustrate  that 
dispersion  is  usually  very  large. 

Other  sections  contain  distribution  information  which  should  also  be 
examined  when  variance  is  considered.  These  other  tables  and  figures  are; 

(1)  The  Square  Law;  Figure  I I — l  (gives  5th  and  95th  percentiles),  II-2,  II- 
3,  (all  the  preceding  are  terminaL  distributions’),  III-L,  1 1 1—2  (the  preceding 
two  are  distributions  functions  of  time-duration,  1 1 1 — 3  (pdf  or  time- 
duration),  VI-3  through  VI-12,  and  Tables  Il-l,  II-2(a),(b),  (preceding  two 
are  terminal),  lll-l(a),  (c)  (time-duration  variances),  VI-3,  VI-4  and  Vl-5 
(the  last  three  are  joint  distributions  at  various  times). 

(2)  The  Square  Law  with  draws;  Table  111-2  (pdf  and  df's  of  time-duration). 

(3)  The  Linear  Law;  TabLe  11-3  (terminal). 

(4)  The  Weale  Special  Model;  Table  II [-3(a)  (time-duration),  VT-b,  and  VI-7 
(last  two  are  joint  distributions). 

Particular  care  must  be  used  in  interpreting  the  figures  on  state 
distributions  which"  are  discrete  as  opposed  to  time-duration  distributions 
which  are  continuous.  That  is,  they  are  all  probability  mass  functions  even 
though  they  are  variously  depicted  as  histograms  or  as  continous.  OnLy  figure 
V-23  is  correctly  portrayed.  The  histograms  should  he  read  as  their  height  at 
mid-points  of  the  rectangles  and  the  continuous  curves  only  at  integers  on  the 
abscissa. 

Also,  the  reader  should  be  aware  of  the  peculiar  manner  in  which  Figures 
V-12,  V-L3,  and  V-14  have  been  prepared.  These  are  frequency  diagrams  (not 
histograms)  and  are  to  be  read  as  having  values  on  1 v  at  Lhe  right  hand  edge  of 
each  rectangle  and  whose  magnitude  is  the  rectangle  height  on  the  left. 

However  the  x(t)  and  y( t )  values  (the  vertical  dashed  lines)  are  located  at 
the  middle  of  the  corresponding  rectangle.  For  example,  in  Figure  V-l 2(a), 
k ( t )  =  12  is  located  at  11.5,  x(t)  =  3  at  2.5  and  so  forth. 

The  examples  shown  in  this  section  amply  support  the  contention  that 
variance  is  Indeed  very  important. 
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WEAlF.  (1972)  pp. 49,50 
SQUAPl  law 


a0  =  bg  =  10 
af  =  bf  =  0 
a  =  0.05,  8  =  0.025 


T 

■A(t) 

™b(  t) 

V  A( t)] 

V(8(t)3 

P* 

0.0 

10.000000 

10. 00001)0 

0.000000 

0.000000 

1.00000 

1.0 

9.756199 

9.606146 

0.494266 

0.744309 

1.00000 

2.0 

9.531694 

9.024 | 77 

o . 

0.477044 

1  .'20000 

3.1 

9.104396 

•.653439 

\ .  4  5  7  6  7  6 

J.  70058  7 

1.00000 

4.3 

9.096330 

3.093493 

’.  . 92  ’  10  5 

0.916767 

1 .70000 

5.0 

9.9001 ’7 

7.64  161- 

7.40642! 

1 . 124; 3.92 

9  #  00990 

6.0 

9. 714510 

/.’O’  ’91 

’ .  3076  ’7 

1  . 279?,39 

0 . 99980 

7.0 

9.539997 

6 . 1 7 1 9 ' 9 

1 . 2  6  ’  6  6  6 

1  .627975 

2 . 99990 

6.0 

9.176300 

6. ’19! ” 

3  .O^T’O : 

1.724010 

0.99466 

9.0 

3 . 7?'1  ]?l 

o .93421  1 

4  *»  c  ^  7  rj 

7  #  970  3*;^ 

0 . 0335 2 

10. 0 

3.0791  ’9 

5.526713 

4 . 36  3714 

’.  1  15356 

0.75673 

1 1 .0 

7.945991 

3.126137 

S  .  39'  10;? 

’.312631 

0.91  ’59 

13.0 

7.322 76  i 

4.731354 

;  #  Q  ’  ;  /  •'* 

2.611096 

0.34826 

n.o 

7. 7093 34 

t. 343692 

r, 

’.7H-><35 

0.75546 

14.0 

7.605639 

5.960360 

’  in’o  jo 

’ .  7? >06 ! 

0.66894 

15.0 

7.511260 

3.6379.90 

■’.717  761 

3.139! 16 

0.5646Q 

16.0 

7.12616? 

3.209579 

3 . 36.9036 

1.362400 

0.45754 

17.0 

7.350737 

2.940190 

°. 957 149 

5.5954  3’ 

0.35954 

13.0 

7.29431  1 

2.471)52 

9.7  737:37 

’..337606 

0.26.391 

19.0 

7.226994 

2.111607 

10.614967 

1.096293 

0.19073 

P*;  probability  mar,; 

not  i h M rb*?d  > 

fin  ®.-itry. 

Table  V-3 


1 


-142- 


CLARK  (1969)  p. 154 

SQUARE  LAW  (WITH  STOCHASTIC  ACQUISITION) 


af  =  bf  =  0 


b0 

a 

a 

PlflOU 

6 

6 

.001 

.004 

5000. 

.284 

1.325 

6 

6 

.0015 

.004 

5000. 

.605 

1.662 

6 

6 

.002 

.004 

5000. 

.944 

1.872 

6 

6 

.004 

.004 

5000. 

1.816 

1.816 

8 

3 

.001 

.004 

5000. 

.202 

1.528 

8 

8 

.0015 

.004 

5000. 

.542 

1.983 

3 

3 

.002 

.004 

5000. 

.966 

2.298 

3 

3 

.004 

.004 

5000. 

2.242 

2.242 

12 

6 

.001 

.004 

5000. 

2.208 

1.879 

12 

6 

.0015 

.004 

5000. 

3.276 

1.743 

12 

6 

.002 

.004 

5000. 

3.662 

1.389 

12 

6 

.004 

.004 

5000. 

2.866 

.400 

Acquisition  Peobabilities  -  (p 

*  q  =  .85 ) . 

Table  V-4 
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CLARK  (1969)  p. 128 
LINEAR  LAW 


af  =  bf  *  0 


a0 

bo 

a 

a 

a 

6 

6 

.001 

.004 

1500. 

.241 

1.346 

6 

6 

.0015 

.004 

1500. 

.530 

1.642 

6 

6 

.002 

.004 

1500. 

.840 

1.811 

6 

6 

.004 

.004 

2500. 

1.678 

1.678 

a 

8 

.001 

.004 

1500. 

.158 

1.571 

8 

8 

.0015 

.004 

1500. 

.441 

1.954 

8 

8 

.002 

.004 

1500. 

.808 

2.192 

8 

8 

.004 

.004 

2500. 

1.990 

1.990 

12 

6 

.001 

.004 

1500. 

1.283 

1.723 

12 

6 

.0015 

.004 

1500. 

2.271 

1.743 

12 

6 

.002 

.004 

2500. 

2.987 

1.493 

12 

6 

.004 

.004 

1500. 

3.120 

.565 

Table  V-5 
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GAFARIAN  and  ANCKER  (1984)  SQUARE  LAW  p.l-12 


[<"*]  A  /» 


o 


II 

o 

a 

r* 

II 

o 


ce 

o 


C<»»r3  */\ 


Figure  V-2 
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Figure  V- 3 
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JU  I  I 


GAFARIAN  and  ANCKER  (1984)  p.27 

NRLQp.321 

SQUARE  LAW 


I  1  ME 


v  [a(0] 


S 


o 

o 


GAFARIAN  and  ANCKER  (1984)  p.I-17 
SQUARE  LAW 


■? 

-  i 

-I 


\ 


3  00 


o  oa 


«  =2,  b  „  =  1 

o  o 


.(=bf=0 


a  =  5 
P  --  80 


SL 


E. 


C  p  :  c, 


E,  Pa*-9 


016  0. 2« 
i  i  roe; 


o  3 : 


o  <o 


Figure  V-A 
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GAFARI AN  and  ANCKER  (19845  SQUARE  LAW  p.I-18 


(O  03 


GAFARUN  and  ANCKER  (1984)  SQUARE  LAW  p.I-20 


HARTLEY  et  at  (1982)  SQUARE  LAW  Fig  3 

a-S  =  o.s 


(0  =  NO.  OF  REPLICATIONS  WHICH  HAVE  NOT  TERMINATED  EARLIER 


FREQ. 


HARTLEY  et  at  (1982)  SQUARE  LAW  Fig  4 

SIMULATION  1300  REPLICATIONS) 

- y<t)  =  NUMBER  SHOWN 


Figure  V-12 


HARTLEY  et  al  (1982)  SQUARE  LAW  Fig  7 

SIMULATION  (SOO  REPLICATIONS) 


a  =  5  =0.5 


Figure  V— 1 A . 


CRAIG  (1975)  SQUARE  LAW  p.29 


CRAIG  (1975)  SQUARE  LAW  p.31 

CRAIG  (1975)  SQUARE  LAW  p  32 

p(a,b,77.6) 

.  pla.b,  116.5) 


Figure  V-L6. 
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CRAIG  (1975)  SQUARE  LAW  p.33 

P(*,b,  f  55,5) 


Figure  V-17. 
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CRAIG  (197S)  SQUARE  LAW  p.57 


Figure  V-20. 
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n 


Figure  V-21. 


HARTLEY  .  HAGUE S  .  AND  t>£  T  HE -WHI  TE  (19U<;) 


c.ye  d  nJ  Lewis  (l  »?o). 


WEALE  (1976)  pp.69,  70 
SQUARE  LAW 


ao  ~  b0  =  10 

af  =  2,  =  3 

a  =  0-5,  B  =  0.25 

P(A)  =  0.776303,  P(B)  =  0.030909,  P(D)  =  0.192788 


DRAWS  OCCUR  IF: 

a  =  2,  b  £  6  or 
b  =  3,  a  <  5 

pfa.b,®)  x  1,000,000 


A  VICTORY 

B  VICTORY 

m 

IKS 

b 

P 

3 

10 

2 

10 

1,506 

3 

9 

199,438 

2 

9 

5,092 

3 

8 

194,456 

2 

8 

9,882 

3 

151,938 

2 

7 

14,429 

3 

DRAM  LEVEL 

DRAW  LEVEL 

■ 

II 

P 

a 

b 

P 

3 

5 

69,044 

2 

6 

17,551 

3 

4 

43,447 

2 

5 

18,646 

3 

3 

26,459 

J 

2 

4 

17,640 

Table  V-9 
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SCALC  OF  OftOINAteS  UNSFFCJFIf 0 

****<1975*1  LINEAR  LAW  p.26 


Figure  V-22. 


WEALE  (1975)  pp.  64,65 
SPECIAL  MODEL 

ao  =  b0  =  10 

af  =  2,  =  3 

P(A)  =  .992809 
P{8)  =  .000018 
P(0)  =  .007173 

DRAWS  OCCUR  IF: 

a  =  2,  b  <  6 
or 

b  =  3,  a  1  5 

ATTRITION  COEFFICIENTS: 

A  SIOE  *  a(0.05  ♦  0.05b) 
B  SIDE  *  b(0.025  «•  .005a) 


p(a,b,<”)  x  1,000,000 


A  VICTORY 

1 

B  VICTORY 

b 

a 

P 

a 

b 

p 

3 

10 

425,276 

2 

10 

0 

3 

9 

336,501 

9 

1 

3 

3 

156,656 

') 

3 

5 

3 

7 

56,612 

2 

7 

i2 

3 

6 

17,764 

L _ I  _ 

DRAW  LEVEL 

DRAW  LEVEL 

b 

a 

P 

a 

b 

i 

P 

3 

5 

5,156 

2 

mm 

27 

3 

4 

1,446 

2 

52 

3 

3 

406 

2 

H 

87 

Table  V-10 


SPRIfWAll  ( 196ft)  p.  170 
SPRIWGACt  MODEL 


1  10,  bg  ■  10,  •  5,  a  *  0.9,  8s  1,  y"  2,  6*2 

tf  a  bf  •  0 


b 

p(o.b.-> 

i 

0.09486 

0.09809 

i 

0.09547 

4 

J.086J9 

* 

*1.0/147 

5 

J.OV/9 

0.03J60 

3 

J.IH/J7 

9 

7.10653 

to 

•1.00135 

SPMNfiAU  ( 19&U)  0.132 

spbiNCAiL  <iooa 


«o  •  15.  bo  •  20,  «2  •  5,  n,  •  5,  a  ■  0.90,  g 


if  *  b^  a  0 


I. 00,  f  •  2.00.  4  •  2.00 


P<«> 

•  0.85503 

P(A| 

*  0.14497 

b 

p(«.b.-l 

4 

P(  4,0,.)  j 

i 

U. 04072 

i 

J  .  7  3*7  '.t i  ; 

2 

0.04855 

i  . . . : » « 1 

) 

0.05657 

) 

>.  124  il 

4 

0.06424 

1 

>.  ll.‘  »6 

5 

J.  U/09b 

•> 

n  *. 

6 

J.  0/602 

•t 

s.M.*U9 

/ 

0. 0/8/2 

9 

7.07850 

1 

i..jo;’9/  i 

9 

7.07504 

9 

1.  1171:  7  1 

10 

0.06838 

:u 

*.«wu/: , 

11 

0.059Q1 

•1.00029  i 

12 

0.04783 

12 

■1. <70010 

13 

0.03602 

13 

O.UOOOJ 

14 

0.02487 

14 

0.00000 

15 

0.01545 

15 

0.00000 

16 

0.00842 

17 

0.00388 

18 

0.00142 

19 

0.00037 

20 

0.00005 

Table  V-ll 
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p(b, 


SPRINGALL  (1968)  MODEL  p.172 
20,  b0  =  20,  nij  =10,  fi  =  3,  QsO.9,  y=2.0,  4=2.0 


1 3  1 4 


1*  20 


SPRINGALL  (1968)  MODEL  p.173 


5,  0  =  5.  m  =  3.  &:30,  0=0.9, 
o  2 


.  4=2.0 


2  3  4  5 


Figure  V-23. 
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VI, 


MISCELLANEOUS  RESULTS 


Two  kinds  of  results  are  presented  in  this  section.  They  are:  (1) 
examples  of  correlation  between  A  and  B  survivors,  and  (2)  some  new  results 
on  combining  small  fire-fights  to  describe  large  battles. 

In  Figures  VI-1  and  VI-2  and  Tables  VI-1  through  VI-7,  the  figures 
provide  visual  evidence  of  correlation  (or  equivalently  covariance)  and  the 
tables  give  actual  computed  values.  The  significance  of  this  is  that 
covariance  (correlation)  is  the  term  whose  existence  (and  magnitude)  causes 
the  Linear  Law  and  any  model  in  which  it  is  a  component  (e.g.  the  Weale 
Special  Model  and  the  Springall  Model)  to  have  differing  L  and  SL  mean  value 
functions.  Other  consequences  of  correlation  have  not  been  explored. 

In  the  remaining  part,  (figures  VI-3  through  VI-12),  two  situations  are 
examined.  Both  include  comparisions  of  L  with  SL  fire-fights  in  which  a  large 
number  of  combatants  either  engage  in  a  standard  Square  Law  battle  or  the  same 
total  initial  numbers  engage  in  several  identical,  simultaneous  and 
independent  Square  Law  fire-fights.  The  two  situations  are:  (1)  48  on  each 

side  for  the  large  battle  and  8  smaller  battles  with  6  on  each  side  and;  (2) 

96  on  each  side  for  the  large  battle  and  16  smaller  fire-fights  with  6  on  each 
side.  The  principal  results  are  that  the  L  model  has  identical  outcomes  no 
matter  how  the  battle  is  sub-divided,  but  the  large  SL  battles  differ 
considerably  from  the  several  smaller  battles.  This  is  true  both  for  the  mean 
value  functions  and  the  standard  deviations.  This  substantiates  the  position 
that  the  compartmentalized  nature  of  real  combat  cannot  be  safely  ignored. 

This  is  some  very  early  work  on  the  problem.  Much  more  needs  to  be  done 
before  any  clear  picture  of  the  seriousness  of  this  situation  will  emerge. 
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WEALE  (1972)  pp.49,50 
SQUARE  LAM 


a0  =  b0  =  10 


af  =  bf  =  0 
a  =  0.05,  6  =  0.025 
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Table  VI-3 
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KARR  (1975b)  LINEAR  LAW  p.12  KARR  ( 1975b)  LINEAR  LAW  p.13 


WEAL E  ( 1976) 


Table  VI-6. 
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SPCCIAt  HOWL 


ANCKER  and  GAFARIAN  (  *  ) 
SQUARE  LAW 


Figure  VI-3. 


(b) 
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ANCXER  and  GAFARIAM  Ofi) 
SQUARE  LAW 


<*A 


-19  9  *•  9  SL 


9a  >•  >*  5L- 


1 9  a  *•  9  ~j  y 
or  SL-^ 

9a  »■  99  j 


SO  100  ISO  200  2S0  300  330  400  4S0  SOO 


ANCKER  and  GAFARIAN  (*) 
SQUARE  LAW 


»o 


p  *p  «.ao 

A  B 


■16  6  V«  a 


SO  too  ISO  200  2S0  300  350  400  450  SOO 


Figure  VI-4 
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ANCKER  and  GAFARIAN  (*)  SQUARE  LAW 


•  *B  3 4 a 

a  9 
Of  sftf  a<J 

yAa«00 


a3.#«  @i  =  * 
33.45  @  Cos 

30.15  <9  1  =  oo 


(a) 


— , - - - r - 1 - 1  —  -| - t - 1 - — ^ 

lOO  150  200  250  300  350  400  450  500 

TIM6 


ANCKER  and  GAFARIAN  (#)  SQUARE  LAW 


5.04 (9t  a  OO 
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THE  VALIDITY  OF  ASSUMPTIONS 
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LANCHESTER  ATTRITION  RATES 
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GIST 


THE  REASON  FOR  PERFORMING  THE  STUDY  was  to  show  that  many  of  the  prevailing 
understandings  concerning  the  relationships  among  the  Lanchester,  Stochastic 
Lanchester,  and  the  General  Renewal  models  of  combat  are  erroneous,  and  to 
collect,  organize,  and  reduce  to  common  notation  almost  all  known  tabled  results 
and  curves  of  particular  examples. 

THE  PRINCIPAL  RESULTS  of  this  study  were: 

(1)  All  Lanchester  model  and  Stochastic  Lanchester  model  mean  value 
equivalent  pairs  differ  for  all  times  except  at  possible  crossing  points. 
These  differences  may  be  very  large. 

(2)  At  least  for  the  Square  Law,  the  Lanchester  model  trajectories  are 
neither  a  universal  upper  or  lower  bound  of  the  Stochastic  Lanchester  model 
mean  value  trajectories. 

(3)  Even  near  time  zero,  the  Lanchester  and  Stochastic  Lanchester  model 
mean  value  trajectories  may  differ  considerably  (they  do  not  differ  materially 
for  the  Square  Law  and  the  sequence  of  one-on-one  duels  version  of  the  Linear 
Law). 


(4)  For  the  Linear  Law,  the  Square  Law,  the  Mixed  Law  and  the  Square 
Law  with  continuous  reinforcements  there  is  a  Law  of  Large  Numbers  on  suitably 
transformed  spaces.  However,  for  untransformed  spaces  this  does  not  apply, 
for  it  can  be  shown  that  as  the  initial  force  sizes  tend  to  infinity  the 
differences  between  Lanchester  model  and  Stochastic  Lanchester  model  mean 
value  trajectories  tend  to  zero,  or  they  may  even  tend  to  a  constant  or 
infinity. 

(5)  Blackwell's  Theorum  does  not  imply  that  individual  combatants  with 
general  interfiring  times  tend  to  have  negative  exponentially  distributed 
interfiring  times.  This  is  even  more  strongly  the  case  for  terminating 
processes. 

(6)  The  Palm-Khintchine  Theorem  does  not  imply  that  superposing  a  large 
number  of  combatants  with  general  interfiring  times  will  yield  a  process  with 
negative  exponentially  distributed  interfiring  times.  This  can  only  be 
approximately  correct  for  large  numbers  and  for  very  large  interfiring  time 
means.  Again,  the  Theorem  is  only  valid  for  non-terminating  processes. 

(7)  Nonhomogeneous  Poisson  processes  do  not,  in  general,  approximate 
general  renewal  processes. 
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(8)  The  Stochastic  Lanchester  model  process  variances  are  generally 
quite  significant  and  can  be  important  for  large  force  sizes,  even  near  time 
zero.  In  addition,  general  renewal  model  process  variances  are  significantly 
different  that  Stochastic  Lanchester  model  process  variances. 

(9)  The  other  Lanchester  model  measures,  (a)  expected  number  of  survivors, 

(b)  expected  time  duration  of  the  battle,  and  (c)  probability  of  winning  are 
even  less  reliable  predictors  than  the  mean  value  trace. 

(10)  The  basic  assumptions  of  the  Stochastic  Lanchester  models,  as  well 
as  the  general  renewal  model,  cannot  hold  for  large  numbers  of  combatants. 

THE  MAIN  ASSUMPTIONS  were  that: 

(a)  All  pre-combat  decisions  have  been  made  and  the  battle  goes  forward 
until  terminated  by  the  action  itself  or  by  tactical  decisions. 

(b)  The  true  model  which  is  to  be  approximated  is  termed  a  General  Renewal 
model . 

(c)  Each  marksman  fires  until  he  is  killed  or  makes  a  kill. 

(d)  The  ammunition  supply  is  unlimited. 

(e)  All  fire  independently. 

THE  MAJOR  RESTRICTION  is  that  the  general  renewal  model  is  a  superposition 
of  many  terminating  renewal  processes. 

THE  SCOPE  OF  THE  STUDY  is  to  show  that  the  basic  assumptions  of  current  combat 
models  do  not  hold  for  large  numbers  of  combatants. 

THE  STUDY  OBJECTIVES  are  to  show  that  many  factors  cause  large  scale  battles 
to  be  a  number  of  simultaneous  and/or  sequential  smaller  scale  engagements 
and  that  use  of  current  deterministic  Lanchester  models  is  incorrect. 

THE  BASIC  APPROACH  is  two  fold.  First  the  Lanchester  model  applications  are 
examined  to  show  fallacies.  Second  all  known  tabled  results  and  curves  of 
particular  examples  are  provided  to  support  the  theoretical  discussion. 

THE  STUDY  SPONSOR  is  TRAC-WSMR. 

THE  STUDY  PROPONENT  is  TRAC-WSMR. 


THE  ANALYSIS  AGENCY  is  TRAC-WSMR. 


